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ABSTRACT 


The alpha=particle decay of states in eae with special 
emphasis on those at 11.83 and 12.71 MeV excitation, has been 
studied using the p+ ¢q nyct? reaction. Neutron time-of- 
flight spectra and solid state detector charged=-particle 
spectra were analyzed by a two-parameter pulse height 
analyzer, Angular distributions and spectrum shapes cof the 
alpha particles from the 12,71 and 11,83 MeV states indicate 
that these states decay sequentially through Be°, The :" 
spin-parity assignment for the 12.71 MeV state is supported 
by the fact that no transitions to the ground state of Be® 
were observed, Approximately 90 percent of the transitions 
from the 11.83 MeV state proceed via Be°(2.9) while the 
remaining 10 percent feed the ground state. There is some 
evidence to indicate that the 10.84 MeV state of eo decays 
exclusively to Be°(0), The 3. assignment for the 9.64 MeV 
state in cl? is favored over 1 on the basis of the angular 
distribution of the alpha particles from this state to the Be® 
ground state, The spectrum shape for this group confirms 
this well established transition. 

This work, together with other work reported in the 
literature, justifies the representation of e as a Be =a 


two-body cluster from the alpha particle threshold to 16 


MeV excitation, Such a representation is in accordance with 
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(Abstract continued) 


the Two-Body Cluster Model of Phillips and Tombrello. 

Appendix II describes a frequency-modulated proton magnetic 
resonance system for the measurement and control of magnetic fields 
from 1 to 12 kilogauss. This system was designed for use with 
the combined broad range spectrograph and spectrometer facility 


of the University of Alberta Nuclear Research Center. 


= F ses 


xt 
/ be: yds is} 
"i ea? a = * 

olferdaeT bee aqtif jdt to (aboM resent? ebe “O1 ide 


alsengea 0 303g bs3s[ubom-yonsupet? s vodixonsh &3 “ breans al 

: : - 

ebleli olssrgen to fortnos bas tnsmeTuasoM etd oT metege os ona 
dghw sey 10? beagiush ssw medaye eid? -aatagor a St 03 & ro} 


2 WY 7 y 
vttitos? 2tsjoemortooqa base dqzrgotioseqe Syria pnotd bsmtdmec etd 


tataeD doinses? «teefoul) adzediA Io yakerev2 


ACKNOWLEDGMENTS 


I should like to express my gratitude to Doctors 
W.K. Dawson, G.C. Neilson, and J.T. Sample for suggesting this 
project, for many discussions pertaining thereto, and for 
much encouragement during the course of the research, 

Special thanks are also due to T.K. Alexander for many 
valuable discussions pertaining to both the research and the 
electronics described herein, 

I should also like to thank Messrs. E. Cairns, J.B. Elliott, 
C. Green and L, Holm, the technical staff of the University of 
Alberta Nuclear Research Center, for much help during the 
course of this work. 

I am very grateful to my wife not only for typing this 
thesis but also for her encouragement during the past three 
years, 

Finally, I thank the National Research Council for 
financial support in the form of three studentships awarded to 


me during the course of this work. 


opm x0 rsbmaxelA .A,T 03 5 ob cae on ih Sanat 
ort bre dovesest of2 daod o2 etn wn ae 
,7ORILA .a,l ,amtted .4 .exeesM aneds of oll opis bluods 7 

Yo ysterevind srt to Msge Ieotaioe2 of alo «J bas aeexD 4D 
of? goiwwh gled dom 36% 283890 doaoeaM xeoloull atzediA 
ie a a 

atd? gmiqys 10? ylmo tom otiw ym of Iulesexg prev me I 
se1ds tesq edd snisub Ineriogetvoons <a | 


a 


o3 bsbraws aeqtdeutnsbuie satis Yo mot ois at 


CHAPTER 


CHAPTER 


APPENDIX 


TABLE OF CONTENTS 


INTRODUCTION 
Motivation 


ara of the Literature on the Disintegration 
of cl 


EXPERIMENTAL APPARATUS 
Introduction 

General Experimental Arrangement 
The Target Chamber 


The Preparation of Self-Supporting Boron 
Targets 


The Fast Neutron Time-of-Flight Spectrometer 
The Charged-Particle Spectrometer 
EXPERIMENTAL MEASUREMENTS 

Introduction 

Pulsed Beam Neutron Time-of-Flight Spectra 
Measurements Using One Solid State Detector 
Measurements Using Two Solid State Detectors 
Conclusions 

THE KINEMATICS OF NUCLEAR REACTIONS 
Introduction 

Programs 

A. Nuclear Reactions I 

B. Nuclear Reactions II 


C. Nuclear Reactions III 


v- 
ton 


eens: 


sR 2 RF RB BB 


wood gnivroqqua-tis® to ona ihe « Pe 
y 


ratamsseq@ sighiTto~seiT noxsusl 3e8T od ww 
reismossoeq@ ototsret-bagnad od a 
aqMsmanueamM JATMDETIESHS 

“ao ksauborsal 

Briosqe right %-YomsmtT sons aod  beetet 
seleuees eevee. bitéu ene ate’ <email 
810399290 94632 bilo? owl anata’ o-teanearea 


Table of Contents (continued) 


APPENDIX II. 


I. 


It, 


IIl. 


IV. 


A TRANSISTORIZED, FREQUENCY MODULATED PROTON 
MAGNETIC RESONANCE SYSTEM FOR THE MEASURE- 
MENT AND CONTROL OF MAGNETIC FIELDS 


Introduction 


Background and Design Considerations 


A General Description of the System 


Circuits 


A. 


B. 


C. 


The RF Oscillator 

The Probes 

The Preamplifier 

The Twin Tee Oscillator 

The Twin Tee Amplifier 

Modulation Amplitude Control Amplifier 
200 cps Compensation Circuit 


Vertical Amplifier for Indicator 
Oscilloscope 


Horizontal Amplifier for Indicator 
Oscilloscope 


Phase Sensitive Detector (or Chopper) 
Driver 


Phase Sensitive Detector 
Time Constant Circuit 


The Power Input Circuit 


Operation 


A. 


B, 


Cc. 


Lining-up Procedures 
Magnetic Field Measurements 


Magnetic Field Control 


98 


98 
99 
108 
ie 
111 
113 
114 
by Be 
116 
117 
| 


118 


118 


Lig 


119 
121 
121 
122 
122 
124 


126 


SS RBRE B 


BI! 


} 


mesey2 oft Se nensetenged: tanenilinll oo 


wosel{tosO TH efT 4A 


1-1 


FIGURES 


Energy level diagram of ct (45 62). Solid 
arrows indicate reported transitions; broken 


arrows indicate unestablished transitions. 
Accelerator drift tube arrangement. 

Target chamber (a) plan view, (b) elevation. 
Electron bombardment apparatus. 

Design details of electron bombardment apparatus, 


Electronics used in the neutron time-of-flight 
spectrometer and in the charged-particle 


spectrometer, 


A pulsed=beam neutron time-of-flight spectrum 
® 
of neutrons produced in the pt ¢a.nycl? reaction: 
BE.” 1 MeV, O = 25°, S = 2,8 meters, Target 
> 0 n 

= 1,2 mg/cm’ isotopic Bll, 

A pulsed=beam neutron time-of-flight spectrum 
produced in the sate nee reaction: E, = 1 MeV, 
07 =o 252s s. = 2,8 meters. Target = 15 agven® 


self-supporting, natural boron. 


Gamma~ray yield curves showing generating volt- 
meter calibration and self-supporting, boron 
target thickness measurement (see text). 


(courtesy of T.K. Alexander) 


Detector geometries for experiments employing 
(a) one solid state detector, (b) two solid state 


detectors. 


Charged=particle-gated neutron time-of-flight 


spectrum showing neutron groups leading to the 4.43, 


7.65, 9.64, 10.84, 11.83, and 12.71 MeV levels of 
12 
Cs 
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07 = 45°, s = 1,15 meters, y = 90°, Ej = 1,5 MeV. 
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(Figures continued) 
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3-10 


3-11 


3-12 


Neutron time-of-flight spectra taken with 
the timing detector in and out of view of 
the cl2 (4.43) recoiling nuclei, 


Charged-particle-gated neutron time-of-flight 
spectrum taken with 10 yin. nickel foil in 
front of the 'timing' detector. oF = 45°, 


s = 1,42 meters, y = 90°, Ey = 1.5 MeV. 


(a) Neutron-gated and direct charged-particle 


spectra 


(b) Charged-particle-gated neutron time-of- 


flight spectrum, e = 45°, s = 1,15 meters, 


y = 90°, E, = 1,5 MeV. 


d 
(a) Neutron-gated and direct charged-particle 


spectra 


(b) Neutron time-of-flight spectrum, OF = 45°, 


sh = 1,15 meters, y = 130°, E, = 1.5 MeV. 


d 


(a) Neutron-gated and direct charged-particle 


spectra 


(b) Neutron time-of-flight spectrum. eT = 45°, 


s = 1,15 meters, y = 240°, E, = 1.5 MeV. 


d 


(a) Neutron-gated and direct charged-particle 


spectra 


(b) Neutron time-of-flight spectrum, ae = 45°, 


S, 7 1.15 meters, y = 265°, E, = 1.5 MeV. 


d 


(a) Angular distributions of alpha particles 
from ct (12.71) to Be’ (2.9) 


(b) Experimental angular distribution of alpha 


particles from co (11.83) to Be (2.9), 
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(Figures continued) 
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3-14 


3-15 


3-16 


3-17 


3-18 


g=19 


Angular distributions for alpha particles 
from ci (9,64) to Be°(0), 


Isometric display of two-parameter data obtained 
for 8. = 45°, i 1,15 meters, y = 90°, 


E = 1.9 MeV, 


Expanded drawing of charged=particle spectrum 
no. 3 of Figure 3-14 with three-body breakup 
distributions. 


(a) Neutron=-gated and direct charged-particle 
spectra 
(b) Neutron time-of-flight spectrum. 9 = 0°, 
n 


sa = 1.15 meters, y = 140°, E.* 1.5 Mev, 


Isometric display of two-parameter data obtained 
for a, = 0°, sh = 1,15 meters, y = 140°, 


7 1.5 MeV. 


Expanded drawing of charged-particle spectrum no, 5 
of Figure 3-17, together with three-body breakup 
distributions, 


Isometric display of two-parameter data obtained for 


@ = 0°, S =1meter, 9 = 90°, 8 = 90°, 
n n Op On 
Ey = 1,5 Mev. 


Isometric display of two-parameter data obtained for 


@ = 45°, S = 1 meter, 9 = 90°, @ = 85°, 
n n On On 


Appendix I 


1 


Kinematics of the reaction pt a nyc? (12.71) for 


E 4 = 1.5 MeV. 


Velocity diagrams for Nuclear Reactions I 


Velocity diagrams for Nuclear Reactions I 
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(Figures continued) 


4 Velocity diagrams for Nuclear Reactions I 

5 Wave vector diagram 

6 Velocity diagram for Nuclear Reactions II 

7 Nuclear Reactions II output for the alpha- 
particle decay of the 9,64, 11.83 and 12.71 
MeV states of ae to the ground and 2.9 MeV 
states of Be° 

8 Velocity diagram for Nuclear Reactions III 

9 Typical Nuclear Reactions III output for the 

11 12% 12% 
reaction B (d,a)C”” «=6((12,71)3; C (12,71) + 
Rx 

ay) + Be’ (2.9) + a, Cy) + a, (8,) + 05(B,)- 

Appendix II 

1 Generalized block diagram of transistorized 
nuclear magnetic resonance system 

2 (a) Rf oscillator 
(b) Preamplifier 

a Probes 

4 Twin-tee 200 cps oscillator 

be Twin-tee amplifier 

6 Modulation amplitude-central amplifier 

7 200 cps compensation circuit 

8 Vertical amplifier for indicator oscilloscope 

9 Horizontal amplifier for indicator oscilloscope 

10 Phase-sensitive detector driver 

BE Phase-sensitive detector 

12 Time constant circuit 
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Figures (continued) 
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16 


Power input circuit 
(a) Front view of NMR control panel 


(b) Back view of NMR control system showing 
positions of the various circuits 


(a) NMR head and probe 
(b) Rf oscillator and preamplifier modules 


Typical proton signals obtained at one kilogauss 
for modulation frequencies of (a) 200 cps; (b) 30 


cps. 


Following 
Page 


127 


ay 


: ie ’ : ne. 
(bauntsqes) 


=e 
~T 


Stusvts duqel rowoT CS 


feneq fo¥anes AMA tq waty ynort = (n) 


qriwods ueseve Londres Mu Yo wetv ase ¢d) 
uiiuotts aeettee eds to anotsteeq 


edoxy ban bead AMM (6) 


Fad 


selybom zefitiqmmetq bus votsitiose a2 «((d) 


aeuagoild sao ts bentaido alsngle aosorg IastqyT 
Of (4) j;uqo OOS (a) 20 eatocerpet? tottaivbom 362 


CHAPTER I 


INTRODUCTION 


I, Motivation 

It is well known that clusters of nucleons such as deuterons, 
tritons, helium-3, alpha particles, and other nuclei in their 
ground and excited states play an important role in the structure 
of light and medium weight nuclei (Wi 58, Sh 60). In particular, 
there is considerable evidence to support the Two-Body Cluster 
Model of Phillips and Tombrello (Ph 60, To 60) in which the wave 
function of a particular state is expanded in terms of a set of 
two-body cluster wave functions, the possibility of three or 
more=body clusters taking part being neglected, This set repre- 
sents all the possible two-body arrangements, A, of the clusters 
that can be formed by the nucleons of a given nucleus, and is 
considered complete and orthogonal for cluster separations greater 
than the range of the strong interactions between clusters. For 
separations less than this range, this set is far from being 
complete and hence the two-body cluster model is not valid. The 
expansion coefficients, or the fractional parentage coefficients 
(62 43 where Rha» gli = 1 under the two-body assumption), are 


energy dependent and are directly related to the probability 
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that the state will decay by the emission of the corresponding 
cluster or, in other words, they are related to the decay widths, 
Explicitly, 

2 2 2 2 
where y is the reduced width and the denominator is the Wigner 
limit (Wi 47). uy is the reduced mass of the two-body cluster, 


and a, is the range of the strong inter-cluster forces. The 


rN 
reduced width is related to the observed decay width ry in the 
usual way by the penetrability factor Py» 
2 
‘ie r,/2P, 
Thus, if for a given two-body channel \', the reduced width eS 
of a state is near the Wigner limit (i.e., if |6 le = 1), it 


2A 
may be inferred that the state is well described by the \' two- 


body, cluster, The measurement of decay widths for particles 

emitted from excited states of nuclei should therefore lead to 

information regarding the two-body cluster parentage of the state, 
Whether or not the whi nucleus fits into the framework of 

the Two-Body Cluster Model is uncertain. Many of the states 


predicted for o™ 


by the more general Cluster Model of Wildermuth 
and Kallenopoulus (Wi 58) using a fe representation and 


average oscillator potentials (Sh 60) have not been observed, 
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However, many of the states that have been observed do fit this 
model quite well. Knowledge of the alpha-particle decay modes 
of the states in o. as implied from the energy spectra of the 


* 
emitted alpha particles, would indicate whether a Be® (or Bee ) 


cluster exists within af. For, if a given level decays predo- 
minantly to the states of Be’, then £18) peSa|? ~ 1 (where s 
denotes the states of Be® involved), and ake state of ci may 
be considered as being primarily of fe ad parentage. On the 
other hand, if the level disintegrates by a direct three-body 
process, the ae two-body cluster is not a good representation 
of the state, 

It is well established that the 7.65 and the 9.64 MeV levels 
of cl? (see Figure l-l for an energy level diagram of oy emit 
alpha particles exclusively to the ground state of aan and that 
the 16,11 MeV level, when it decays by particle emission, pri- 
marily feeds the 2,9 MeV level of Be® although transitions to 
the Be® ground state and some three-body emission have also been 
observed, The levels between 13 and 16 MeV excitation have been 
observed to decay predominantly to Be°(2.9), It would seem then 
that these levels do fit into the two-cluster picture of a 
and indeed the spins and parities of the 7.65, 9.64 and the 16.11 
MeV states can be explained in terms of the cluster model of 
Wildermuth and Kallenopoulus (loc, cit.) using the —— repre- 


sentation (see reference Sh 60). 


Very little or no work has been reported on the alpha- 
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Figure _l-1, Energy level diagram 
of Cl2 (Aj 62). Solid arrows 
indicate reported transitions; 
broken arrows indicate unestab- 
lished transitions, 
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particle decay of the levels between 9.64 and 13 MeV; namely, the 
10,1, 10.84, 11.83, and the 12.71 MeV levels, Since there is 
speculation (see Chapter III) that the 12.71 and the 11,83 MeV 
levels might decay by a direct three-body process, in contrast to 
the sequential process through ne® predicted by the Two=Body 
Cluster Model, an experimental program was initiated to inves= 
tigate the alpha=-particle decay mechanism of these states, That 


program constitutes the subject matter of this thesis, 


Lia Survey of the Literature on the Disintegration of cl? 


The disintegration of er has been the subject of experimen- 
tal and theoretical research since 1934 when Chadwick, Feather, 
and Davies (Ch 34) observed, in a cloud chamber containing methane, 
a three-pronged star among tracks of protons projected forward by 
radium~beryllium neutrons, They attributed this event to the 
reaction, brane. 

Dee and Gilbert (De 36) in 1936 provided the first successful 
analysis of a nuclear reaction in which more than two reaction 
products were involved, explaining the energy distribution of the 
alpha particles from the reaction 5 | (o-wie, In this case, rr 
at about 16 MeV excitation was involved, the proton bombarding 
energy being 200 keV. The alpha-particle spectrum was divided 
into three groups: (a) a homogeneous group at 5.7 MeV, containing 
about one percent of the alpha particles, and attributed to the 


Bo Ya oN ite weod) reaction; (b) a broad group at 3,85 MeV, of 
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width 0.77 MeV, attributed to the Bn tas iper aus MeV) reaction; 
(c) a continuous group extending from low energies to about 5 
MeV, containing twice as many particles as group (b) and attri- 
buted to the alpha-particle breakup of the first excited state of 
Be®, 

In 1949, Hanni, Telegdi, and Zunti (Ha 48) discovered the 
photodisintegration of carbon in nuclear emulsions when they 
observed three~pronged stars in plates which had been exposed 
to gamma rays obtained from the bombardment of lithium with 
protons, The analysis of these stars led to the result, similar 
to that of Dee and Gilbert (De 36), that the breakup of haa 
proceeded predominantly through the first excited state of i: 
There was no evidence in this work for any direct three-body 
decay of eo", 

Since this work, many experimenters have investigated the 
oa Oe reaction (Te 50, Go 50, Te 51, Wi 53, El 52, Go 53, 
Li 53, Mi 53, So 55, Ha 55, Ma 58, To 60a, Se 61) using nuclear 
emulsion techniques at gamma-ray energies up to 300 MeV. The 
cross section shows two broad peaks at gamma~ray energies of 
19 and 29 MeV, and a pronounced minimum at 21 MeV, the region 
of the (y,n) giant resonance, It is zero above 42 MeV, For 
gamma-ray energies below 19 MeV, the decay of ee? proceeds 
predominantly through the first excited state of Be°(2.9 MeV), 


and for energies above 26 MeV the reaction goes via the ground 


8 
state of Be® (12 percent) or through the states in Be in the 
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region of 17 MeV excitation, In the region between 19 and 26 MeV 
in gamma-ray energy, the number of events which go through the 
ground state of her increases with energy from 5 to 18 percent, 
Decay to the 2.9 MeV state can account for at most one-third of the 
events in this region. It is uncertain whether higher states in 
Ba® or direct three-body decays are involved in the remainder of 
the events. 

Many other reactions and experimental techniques have been 
employed to investigate the mechanism of the breakup of gi into 
three alpha particles, These will be mentioned and discussed in 
the following summary of the information available on this 
process for each level of o from the breakup threshold (7,281 


MeV) to the 16.11 level. 
(a) The 7.656 State 


1, Miller, Rasmussen, and Sampson (Mi 55) investigating the 
12 Pe Yh 
Cre (ayatnc reaction with 22 MeV alpha particles failed to 
observe the ci? recoil corresponding to the alpha group leaving 
c’ excited to this level, They estimated that the probability 


of the decay proceeding through Be® as greater than 80 percent. 


2. Bent, Bonner, McCrary, and Ranken (Be 55) using the 
Bes(aences. reaction to feed the 7.65 state and a magnetic lens 
pair spectrometer to study the subsequent radiations, concluded 
from the relative intensities of the 4.4 MeV gamma radiation and 


the internal pair line that the 7.65 state breaks up into 
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Be? +a with a probability greater #han 96 percent. 


3. Cook, Fowler, Lauritsen and Lauritsen (Co 57) were the 
first to observe the alpha-particle decay of the 7,65 level, They 
measured the alpha-particle spectrum associated with the 8 decay 
of Be formed by the deuteron bombardment of Bet and concluded 
that the 7,65 MeV state breaks up into alpha particles via Be” 


and that the results are inconsistent with a direct three-body 


decay process, 
(b) The 9.64 State 


1. Livesey and Smith (Li 53) reexamined photographic plates 
exposed to neutrons from deuterons on lithium by Green and Gibson 
(Gr 49), The energy spectrum of the inelastically scattered 
neutrons showed two peaks corresponding to energy levels in 
at 9,6 and 11.8 MeV, The events which were attributed to the 
9.6 level were consistent with the assumption that they involved 


the ground state of te The events attributed to the 11.8 level 


went through the first excited state at 2.9 MeV. 


2. Jackson and Wanklyn (Ja 53) investigated the disintegra- 
tion of ci using fast neutrons (E. up to 45 MeV) from the proton 
bombardment of fe. and cloud chamber techniques. For ED less 
than 20 MeV, most of the events proceed to the ground state of 


her through a level in pag at 10 * .& Mev. 


3. Frye, Rosen, and Stewart (Fr 55) irradiated C-2 emulsions 


with neutrons of fifteen discrete energies in the range 12-20 
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MeV (obtained from the T(dn)He’ reaction), and found evidence for 
the excitation of the 9.6 level in av and the ground and 2.9 states 
in Be’, Six events appeared to have gone through the 7.65 state in 


ot, 


4, Perkins (Pe 51) also investigated the disintegration of 
ct by neutrons, He irradiated nuclear emulsions with neutrons 
(EO up to 24 MeV) obtained from the Li? + d reaction, All 485 stars 
examined indicated that the disintegration proceeded through fe", 


No particular levels of ge were specified, 


5. Vasil'ev, Komarov, and Popova (Va 58) exposed nuclear 
emulsions to neutrons of energy up to 19 MeV obtained from the 
bid.) reaction. In contrast to other work on the inelastic 
scattering of neutrons from e- (Pe 51, Ja 53, Fr 55), the research 
of these experimenters indicates that a large fraction of the 
disintegrations occur without the formation of any intermediate 
nuclei, For example, in the neutron energy range from 9 to ll 
MeV only 10 percent of the events involve the e nucleus (below 
9.6 MeV), and only 40 percent of these events involve bas These 
values are typical for the entire neutron energy range studied. 
The energy distribution of the reaction products follows a 


statistical (phase space) distribution. 


6, Need (Ne 55) bombarded methane in a cloud chamber with 
12 
29 MeV protons to investigate the mechanism of the C’ (p,p')3a 


reaction, He excited levels in cl? near 9,6, 12, 16, 20, and 25 
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MeV and found that the 9.6 level decayed entirely to the ground 
state of Bay The level in the vicinity of 12 MeV went to either 
the ground state or the first excited state of Be® with about equal 
probability. The levels between 13 and 16 MeV were found to decay 
to the ground state of ne” in less than five percent of the cases, 
with the rest going to the 2.9 state of a, For the levels at 

20 and 25 MeV, it was estimated that 16 * 9 percent and 27 + 10 
percent respectively of the events involved He” in its ground state, 
and that between one-third and one-half of the remaining events 
proceeded through the 2.9 level in hee The mechanism for the 
remaining events was uncertain, There was no evidence of direct 
four=body decay of ae In similar experiments, Vasil'ev, 
Komarov, and Popova (Va 63) using protons in the energy range from 
15 to 29 MeV found that the primary mode of decay is through the 

N = compound system to eg and Ree in the ground, 2.9, and 

11.3 MeV states. Another group of Russian experimenters, Vasil'ev, 
Komarov, Koshelyaev, and Popova (Va 60), observed direct interac~ 
tions of the p-a type in their measurements on the angular dis=- 
tribution of alpha particles from the reaction oul COLD for 


a proton energy of 29 MeV. 


7, LaSalle, Bent, and Cramer (La 63), using two solid state 
detectors and two-parameter pulse~height analysis to look at 
fast coincidence events between inelastically scattered alpha 
particles (bombarding energy of 22.5 MeV) and those resulting 


from the breakup of unbound levels in ig concluded that the 
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disintegration of the 9.64 level goes primarily by a two-body 
process through Be’, LaSalle, Cramer, and Eidson (La 63a) 
measured the angular correlation of the alpha particles from the 
breakup of the 9.64 state to the ground state of oi and con= 


firmed the spin assignment of 3 for this state in ee. 


&, Vasil'ev, Komarov, and Popova (Va 62) also initiated 
the disintegration of c* using alpha particles at 23 MeV, 
They conclude that the decay of c}? into three alpha particles 
proceeds with high probability through the direct decay of the 
compound system o. The intermediate nucleus, te” appears 
as the result of the resonance interaction of the final state 


alpha particles. 


9, Bogatin, Novak, and Ostroumov (Bo 63) used nuclear 
emulsions to study the breakup of grF induced by the inelastic 
scattering of 80 MeV nt mesons, The excitation of the 9.6 MeV 
level was found to contribute about 20 percent to the reaction 
cross section, The energy distribution of the alpha particles 
from this state was explained on the basis of the assumption 
of a resonance interaction between the final state alpha par- 


ticles. 


10, Work on the photodisintegration of oxygen (Da 56, 
Aj 59) indicates that for gamma-ray energies below 22 MeV, 
the reaction takes place mainly through the 9.64 and the 10.8 


8 
MeV states of co to the ground state of Be . 
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11, Brown and Knowles (Br 62) from work on the alpha- 
particle bombardment of cl? using 25 MeV projectiles confirmed 


that the 9.64 MeV state decays exclusively to Be’, 
(c) 10.1 State 
1. See item (b)2. 


2, Cook, Fowler, Lauritsen, and Lauritsen (Co 58) observed 
the formation of the 10.1 state in e- in the 8. decay of hae 
More than 96 percent of the events corresponding to this level 
decayed to the ground state of acs with Be? (2.9 MeV) involved 
in the remaining events. Direct three-body disintegration was 
ruled out on the basis of the shape of the alpha-particle energy 


distribution, 


3, Wilkinson, Alburger, Gallmann, and Donovan (Wi 63) 
more recently confirmed the work of Cook et al (Co 58) in their 


investigations into the alpha=-particle emission following the 


beta decays of pi? and a, 


(d) The 10,84 State 


No work has been reported on the alpha=particle breakup of 


this level other than that cited in item (b)2 and item (b) 10. 
(e) The 11,84 State 
1. 3See itema(b)1. 
(f) The 12.71 State 


1, See item (b) 6, 
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2. Vasil'ev et al (Va 63) found evidence for the formation 
of a level near 12.7 MeV in gt; in their work on the proton bom- 
bardment of carbon, for proton energies in the range 15 to 20 
MeV. There was also evidence for the 0, 2.9, and the 11.8 MeV 


states in Be®, 
(g) States Between 13 and 16 MeV 
iy See item (b)6. 
(h) The 16,11 MeV State 


1. The ete ane (2.9 MeV) is a very well established 
reaction and was in fact one of the first of all nuclear reac- 
tions to be investigated, Dee and Gilbert (De 36) were the first 
to explain the breakup of this level on the basis of the existence 
of an unstable fig nucleus of very short life, as mentioned 


earlier. 


2. Geer, Nelson and Wolicki (Ge 55) assigned spin-parity 
values of gt for both the 16,11 MeV state of gh and 2.9 MeV 
state of Be® on the basis of angular correlations between the 


alpha particle leading to Be? (2.9) and one of the alpha particles 


from the breakup of this level. 


3. Beckner, Jones, and Phillips (Be 61) employed a 
magnetic spectrometer to resolve the reaction products produced 
in the B’(p,a)Be™ reaction, and compared the alpha-particle 


continuum with the shape predicted by the generalized density~- 
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of-states function of Phillips, Griffy, and Biedenharn (Ph 60a) 
which is based on the assumption that the three-body decay 
proceeds sequentially through ge The comparison is good 


and therefore seems to confirm this mode of decay. 


4, Dehnhard, Kamke, and Kramer (De 61) suggest a simul- 
taneous breakup of this level within the resonance a = 163 
keV), and confirm the sequential process through te” for 


proton energies immediately above and below this resonance. 


5. Bronson, Simpson, and Phillips (Br 63) investigated 
this reaction using two solid state detectors and two- 
parameter analysis techniques. Dalitz plots of the energy 
distribution showed peaks corresponding to known levels in 


Be’, in particular, the ground and 2.9 levels. 
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CHAPTER II 
EXPERIMENTAL APPARATUS 


I. Introduction 


The alpha-particle breakup of the 12.71 and the 11.83 MeV 
levels in ed fed by the rete Mees han reaction, was studied using 
neutron time-of-flight techniques and charged particle spectrometry, 
A self-supporting, natural boron target of thickness 15 uemafem™ 
was bombarded with 1.5 MeV deuterons obtained from the University 
of Alberta 2 MeV Van de Graaff accelerator. The states of ree were 
selected by a neutron time-of-flight spectrometer designed by 
Neilson et al (Ne 59), and the alpha particles from the decay of 
these states were detected by solid state surface barrier detectors 
located in the primary reaction plane. Zero-time reference pulses 
for the neutron spectrometer were generated from the detected 
associated alpha particles by means of a fast electronic system 
developed by Alexander and Neilson (Al 63). 

Two experiments were performed in this investigation. In the 
first, one solid state detector was employed, and the charged 
particle and the neutron time-of-flight spectra were analyzed by 
a two-parameter, 1024 channel, TMC* pulse-height analyzer operating 


in the 32 x 32 mode, and a single~parameter 256 channel, TMC* unit. 


* Technical Measurements Corporation, 441 Washington Avenue, 
North Haven, Connecticut 
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In the second experiment, two solid state detectors viewed the target, 
One of these provided the timing pulses while the other, gated by the 
"timing" detector, provided coincidence particle spectra which were 
also analyzed with the neutron time-of-flight spectrum by the two= 


parameter analyzer. 


I, General Experimental Arrangement 


The deuterons for this experiment were obtained from the 
University of Alberta 2 MeV Van de Graaff accelerator manufactured 
by High Voltage Engineering Corporation, This machine provides an 
energy stability of about +3 kev. It is located, together with a 
beam deflecting system and the analyzing magnet, in a room with 24 
inch thick concrete walls to keep gamma and neutron background to a 
minimum, The accelerator drift tube arrangement is shown in Figure 
2-1. To provide additional background reduction, the drift tube is 
encased in a one inch thick lead sheath from the point where it enters 
the target room to approximately one meter from the target chamber. 
Additional lead and 4 inches of a paraffin-lithium carbonate mixture 
surround the carbon trap, the energy stabilizing slits, and a retrac- 
table beam viewer (not shown). Because the target room is of light 
construction and has a fifteen foot square and seven foot deep exca~ 
vation under the floor at the target position, neutron scattering is 


kept to a minimum. 


The neutron detector carriage was mounted on a boom which 


pivoted about a point directly below the center of the target chamber. 
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The neutron flight path was adjusted by sliding the carriage along the 
boom, and the angular position of the detector was read from a protrac=- 
tor. The solid state detectors were mounted inside the target chamber 


which will be described in the following section. 


Lil. “The Target Chamber 


The target chamber used in these experiments is shown in plan 
and elevation views in Figure 2-2, Its design was conceived by the 
chief technician of our laboratory, Mr. J.B. Elliott. 

The target was held in position at the center of the chamber in 
a stainless steel frame A, whose vertical and lateral position as well 
as angular position in the reaction plane (the horizontal plane) 
could be externally adjusted, The mechanism providing this adjustment 
was mounted on the half inch thick teflon base plate of the chamber, 
thus insulating the target from ground, 

A brass Faraday cup, B, mounted on a brass rod which passed through 
the base plate was used to collect the beam after it had passed through 
the target. Actually, the target holder and the Faraday cup were 
connected together so that a true current reading would be obtained. 

Beam alignment on the target was facilitated by the quartz 
viewer, C, also mounted on a brass rod passing through the teflon 
base plate, It could be positioned directly in front of the beam 
entrance to the chamber. 

Provision was made to allow two solid state detectors to view 


the target. One of the detectors was mounted in a transistor socket 
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Figure 2—2. 


Target chamber (a) plan view 


(b) elevation. 
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on the end of the stainless steel tube D, which slips inside the brass 
tube E, An O-ring was used to effect the vacuum seal between these 
two tubes. The brass tube was welded to a stainless steel strap F, 
1.25 in, wide and .010 in, thick, which provided a sliding vacuum 

seal over a .5 in. wide slot cut in the wall of the chamber. This 
arrangement permitted the continuous adjustment of the angular 
position of the detector from zero to #140° with respect to the beam 
direction. The signals from this detector were brought to a BNC 
connector mounted on the end of the stainless steel tube D, 

The other detector was mounted in a holder, H, fastened to but 
insulated from a rotatable brass sleeve J which was part of the 
target-holding assembly. Thus, continuous adjustment of its angular 
position was also available. The signals from this detector were 
brought out of the chamber by means of a brass rod that passed 
through the teflon base plate. Both detectors could be set to 
approximately +1/4 degree, and electrical grounding was made to the 
body of the chamber. 

During the course of an experiment, both the top and the 
bottom of the chamber were electrically shielded by aluminum covers 
In order to facilitate the calibration of the solid state 
detectors at any time during the course of an experiment, a * 
alpha-particle source was mounted inside the chamber on the end of 


another brass rod, G, whose position could be adjusted from outside 


the chamber, 
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IV. The Preparation of Self-Supporting Boron Targets 


(a) Introduction 

The use of solid state counters to detect nuclear reaction 
products necessitates the use of targets which have little or no 
backing material, The elastic scattering of the incident beam into 
the counter by a thick, or thickly backed, target would cause 
serious electronic pile-up and would unduly shorten the lifetime of 
the detector, Thin self-supporting targets are also desirable from 
the point of view of keeping the energy degradation of the reaction 
products to a minimum, thus facilitating the identification of the 
various reactions taking place. 

Self-supporting boron (natural a and separated 2’) targets, 
3/8 in, to 1/2 in. in diameter and ranging from 10 to 25 denen 
were prepared by vacuum evaporation techniques using electron bom- 
bardment heating. 

In electron bombardment heating, a focussed beam of electrons 
impinges directly on the surface of the material to be evaporated, 
This is a very efficient method of vapour source heating, and is 
superior to most other methods because the highest temperature 
region is the surface of the evaporant itself rather than the 
crucible supporting the evaporant. Thus by maintaining the tem- 
perature of the crucible-evaporant interface at a low level, it is 
possible to reduce the degree to which an active evaporant attacks 
the support material, prolonging crucible life, and keeping to 


negligible proportions target contamination by the support material. 
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Figure 2—3. Electron bombardment apparatus 
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Resistance heating of the boron using crucible~filaments made 
of refractory metals such as tantalum, tungsten, or molybdenum proved 
to be impossible because of the reactivity of the boron with these metals 
at the temperatures required for the volitilization of the boron 
(2300°C). The use of carbon crucibles was ruled out because of the 


carbon contamination of the boron targets which would have resulted. 


(b) The Electron Bombardment Apparatus 

The electron bombardment apparatus was a modified version of a 
design given by Muggleton and Howe (Mu 61), and is shown in Figure 
2-3. Figures 2-4a and 2-4b contain the design details of the parts 
involved, The apparatus was mounted on a Véerd | vacuum system, 
model VS-9, which is equipped with a two inch oil diffusion pump 
and a rotary forepump, The pumping speed provided by this system 
is satisfactory for resistance heating evaporation, but is just barely 
adequate for electron bombardment evaporation. Because of this 
fact, the material in the pumping path from the crucible to the 
vacuum system was kept to the absolute minimum. 

With reference to Figure 2-4a, it can be seen that this type 
of heating device is basically a simple diode. It was mounted on 
a 14 in, diameter, 1 in, thick, stainless steel base plate and 
contained in a 12 in, diameter bell jar. The anode was formed by 
the tantalum crucible A and its tungsten support rod B which was 


ts 
insulated from the base plate by the Kovar insulator C. The 


* Vacuum Electronic Engineering Co. 


xk Available from Stupakoff Ceramic & Manufacturing Co., Latrobe, Penn 
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Figure 2-4(a). Assembly drawing of electron bombardment 
apparatus, 
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Support rod was a sliding fit in the brass boss D to provide easy 
height adjustment of the crucible, A glass insulated lead was used 
to connect the anode to one of the type 95.0056 Kovar terminals 
mounted on the base plate, An NJE™ 5 kv, 0-500 milliampere, dc 
power supply was used as the high voltage source for this system, 
The filament (cathode), F, consisted of a 3 in, x 1/16 in. strip 
of ,005 in. thick tantalum, spot welded to two tungsten rods E, 
and shaped as shown in Figure 2-4b with a 1/2 in. diameter loop 
in the center, The two filament support rods were held rigidly 
3/4 in. apart by the ceramic spacer, G, and were connected to two 
of the Kovar terminals on the base plate through the bars H. The 
filament current supply was rated at 12 KVA, or 50 amperes at 24 
volts, 

A water cooled heat shield, I, was necessary to protect the 
bell jar from the intense heat generated during the evaporation 
process, It was also found essential to provide water cooling, J, 
to the substrate onto which the boron was to be evaporated to 
prevent the deposited film from rupturing. 

Electrons emitted from the filament were focussed onto the 
surface of the boron in the crucible by the molybdenum focussing 
electrode K which rested on a ledge machined into the heat shield. 
The size of the aperture in this electrode, the diameter of the 
filament, and the geometry of the anode-focussing electrode- 
filament system determined the electron spot size and the dc 


voltage and current required to achieve the evaporation temperature, 
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(c) Procedure 

A nickel slide, 3 in. x 1 in, x .020 in., cleaned in a concen- 
trated sulphuric acid-chromic acid cleaning solution and rinsed in 
distilled water, was used as the substrate for the evaporated boron. 
To facilitate the subsequent stripping of the film, the slide was 
coated on one side with a thin layer of a liquid detergent, Teepe. 
and polished vigorously with grease free tissues, The treated slide 
was then placed into the water cooled holder, Teepol side down, and 
copper shims were inserted behind the slide to ensure good thermal 
contact with the cooling jacket, All handling of the substrate was 
done with tweezers, 

Prior to the positioning of the substrate, the tantalum 
crucible was cleaned in the same manner as the slide, and filled 
with boron, The diode geometry was then set, Best results were 
obtained with the diameters of the aperture of the focussing elec- 
trode and the filament fixed at 1/2 in. The (top of the) crucible 
and the filament were positioned 1/8 in. below and 1/4 in, above 
the focussing electrode respectively. The spot size corresponding 
to this geometry at the evaporating conditions of voltage and 
current was approximately 1/8 in. in diameter. To ensure reasonable 
uniformity of the boron deposit over the slide, the distance from 
the vapor source to the substrate was permanently fixed at 4 inches, 

The evaporation chamber was then pumped down to a pressure of 


approximately 1 x yore mm of mercury, Before turning on the 


filament current, a 3 in, diameter shadow shield was placed between 


* Teepol Liquid Detergent -- Available from Shell Oil Co. 


a —_——— SS = 


- iS = 


<pSanos 2°ni bennslo ,.aF O80. x sat SD x eat © yebtle® 
dt beets bose coksufoe qainssia bios atmewto-bies ote 
notad begaxoqevs sds rol siatdedoe eds a6 boeu easy . 7838) 
enw obife si7° mi!) oda Yolgnbqghite Idsopoedue 
*pemadt Ikegtageb biupil s to tevst aids 8 diiw obte 
shite borers SdT .Aguent? sext stestg dotw ylevongtv 
bas ,mwob eshte JoqeaT .rebLor baloeo agzew 9nd gaat 6 
[ntted? boon ovveds of abtlz add bnkded betresnt e2pw 
anw stovtiadue, ert to aatfbnsnd EPA -redoat gnifoor add 
 eTOReaN 
mufatnsi od? ,sosttedve sit to aataolsteoq ada oF : 
borll? bre ,oblfe off es 19am sme of3 nt bonnes enw | 
sysweesivasy gest se nod asw ¢xdamosg both aM ynoxed dh Ae 
“5929 antago? oft to stvtseqs sf? Yo etezemsth ia | aie 
aldterts (ed2 Yo oot) oct Vat s\s 38 beet? srome Lid ads 
svods ,nit $\{ boa woled ynt &\E bano2steog o19m sient 
anibnoqastyoo este toga sit svievisoeqeoy shoe 
fam agntioy to. ea rot dtbeas guidsroqsve Sd2 ds 


sfdknoasst atuere ol ,1939melb af -nl B\! ylesamh 
Mott sansjath odd ,4bita aie ave —— 
-ratont & 38 bextt vitenoianhe bay pea te 
io otvarsty 5 of mull beqmuq ‘aap enw ° 
odd a0 onli, stoted : ene \ a 


tpwd od haantry AoW bTabity, atyatlae 


Pr _ 


j a 


109 LY0 I fous mos) 9 


oF wn 


the substrate and the crucible to minimize the deposition of conta- 
minants on and to avoid undue heating of the substrate during the 
outgassing process, 

The filament temperature was then raised to approximately 
1500°C, as determined by an optical pyrometer, by setting the 
filament current to approximately 20 amperes. After allowing about 
15 minutes for the outgassing of the filament, the positive voltage 
was applied to the anode and gradually increased to 500 volts. A 
condition of space charge limited emission from the filament was 
established by adjusting the filament current to a value above which 
there was no further increase in the anode current. These conditions 
were maintained until the boron and the crucible were outgassed 
which, for a fresh boron charge, took the order of twenty minutes. 
The high voltage was then very slowly increased, maintaining space 
charge limited emission and never allowing the pressure rise above 
2x 107" mm of Hg, until the boron started to evaporate. With the 
geometry used, a satisfactory rate of evaporation took place at a 
voltage setting of 3000 volts and an anode current of 100 milliam- 
peres, Overheating was avoided as it caused spitting of the boron 
which resulted in pinholed films and undue loss of boron from the 
crucible, 

When the above evaporation conditions were reached, the nickel 
substrate was exposed to the vapor source, Exposure times for films 
between 10 and 25 Noha feme thick ranged from 5 to 10 minutes. 


Thicker targets could be made by increasing the exposure time, of 
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course; however, films thicker than about 25 ugitefen® were very 
difficult to strip as they tended to curl up when released from the 
substrate, 

The pressure of the system during the evaporation period 
should be kept below 1 x 107° mm of Hg to avoid discharges from 
the anode to ground, 

At the end of the evaporation, the high voltage and the 
filament current were turned off, and the system was allowed to 
cool for about ten minutes before being vented to the atmosphere. 

The boron film was released from the substrate by slowly 
sliding it into a solution of aqua-regia (one part water, one 
part concentrated HCl, and three parts concentrated HNO ,) at a 
very small angle with respect to the surface of the solution, This 
solution etches the nickel at the boron-nickel interface without 
attacking the boron film. Prior to this releasing procedure, the 
film was scored into 1 in. squares, and the edges of the substrate 
were rubbed with a wooden pencil in order to loosen the film at its 
edges. When the first section of the film was released, the 
substrate was removed from the acid solution and dipped into dis= 
tilled water at an angle of about 30 degrees with respect to the 
water surface, The surface tension of the water lifted the film 
cleanly from the substrate. Care was taken to prevent water from 
creeping over the top of the film so that it would not curl up or 
sink when it became free of the substrate. 


Before releasing the other sections of the film, the first 
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one was mounted on a target frame, The frame, which was a 1 1/4 in, 
x 3/4 in. x 1/32 in. piece of brass or nickel plate with a 3/8 in. or 
1/2 in, diameter hole in the center, was placed under the floating 
film and raised vertically, allowing a small strip of the film to 
wrap around the edge of the frame, The edges of the aperture in the 
frame were highly polished to prevent the film from tearing on sharp 
corners, Excess water was carefully blotted from the frame before 
it was placed in a draught-free atmosphere to dry. The finished 
target was then stored in a vertical position in a tightly closed 
plastic box. Mechanically, the targets were extremely fragile when 
wet but became more rugged with drying. However, extreme caution 
was always required when they were handled. 

These targets, because boron is practically a refractory 
material itself, were able to withstand a considerable amount of 
bombardment by the beam of the Van de Graaff accelerator. For 
example, one target withstood a 0.2 yamp beam of 1.5 MeV deuterons 
for the equivalent of about one week's continuous running before 
it had to be replaced, 

In general, the purity of the targets was very good. Carbon 
and oxygen were the chief contaminants, and even these were just 
barely detectable in some targets. The degree of oxygen contamina~ 
tion depended on the length of time the film was exposed to the 
aqua-regia stripping solution, This time was therefore kept to a 
minimum. 


The thickness of the target used in most of the experiments 
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was measured by T.K. Alexander of this laboratory using the 873 keV 
resonance in the Festa yor’ reaction, The gamma-ray yield was 
measured with and without the boron target in front of a thick 
fluorine target as a function of the proton bombarding energy, 
Figure 3-3 shows the two yield curves obtained, The thickness of 
the boron target in keV was taken as the shift in the position of 
the low energy side of the yield curve, which was in this case 5 
keV for the 873 keV protons, This corresponds to approximately 


L5 gts ‘eine of boron, 


V. The Fast Neutron Time-of-Flight Spectrometer 


The neutron time-of-flight spectrometer used in the experiments 
to be described in this thesis was designed by Neilson et al (Ne 
59, Ne 60, Ne 62), Because its "start" pulses were obtained in a 
manner not treated in detail in these references, and for the sake 
of completeness, a brief discussion of the instrument will be given 
here, 

This spectrometer was designed to determine the energy of 
neutrons in the MeV range where the flight times are of the order 
of tens of nanoseconds, To this end, it detects and measures the 
times of flight of neutrons produced in nuclear reactions, the 
flight path being defined by the distance between the target and 
the neutron detector. For each time measurement there are two 


pulses required; a "start" or zero~time reference pulse which is 


linearly related to the time at which the neutron generating reaction 
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takes place, and a "stop" pulse which indicates the arrival of the 
neutron at the detector. A time-to-amplitude converter is used to 
produce an output whose amplitude is proportional to the extent of 
the overlap of these two pulses, In other words, the output is 
proportional to (1 - t/T), where T is the length of the "start" and 
"stop" pulses, and t is the delay between the front edges of these 
pulses and, of course, is a direct measure of the neutron flight-time. 
The neutron energy (nonrelativistic) is simply related to the flight 


time by the formula, i * fe a 


nanoseconds/meter, where EO is 
the neutron energy in MeV. 

There are three methods by which the "start" pulses can be 
derived for this instrument. Two of these involve detecting the 
associated gamma-rays and charged particles resulting from the decay 
of states, bound and unbound respectively, formed in stripping 
reactions. The third method employs beam pulsing techniques and can 
be used independent of the type of reaction producing the neutrons, 
Here the bombarding beam of particles is pulsed onto the target for 
very short periods of time, typically one nanosecond, at a rate of 
1.67 mc/s. The method used to pulse the beam is described in detail 
in the above references but will be mentioned here briefly because 
this technique was employed in some preliminary measurements of 
this work, 

The zero-time reference pulses are derived from a signal induced 
in a pick-up loop by the sinusoidal voltage used in sweeping the dc 


beam of the Van de Graaff accelerator across the "scraper" slits 


~ os * 


edt to faviass oft astsotbrl dott ealuq “qote" a bee, 28. 4 


03 baeu al redtevKn0o abut iqua-os~smkd A ToJoo38b 
" 90 share 


to anes ef9 ad Ssaotixoqotq et ebustiqnua saodv segues 


rr 
at tuqivo ofa ,ebtow vedio al -eealuq ows assis to galzav vo edz 


latin Sate iy 


bas “Sende” oda to dagnel ody et T exonw .(T\3 - i) oF fanots; 2 


Si 
seeds Yo aenbs iaoxt eds meswted yafeb oft at 9 baw — whe : 
4 en . 
ae sotiven sd lo sxuesom toetlb s at ,semu0s to ,b 
gigtli edt of betnles yiqute at (otveatvioslermon) ygrede a 
at 8 erodw ,redom\ebroos20nsa SM ave St 2 4 ,siverto? pe: an + 
.VeM nt vgrens aosauen of 9 


sd ag> eseiuq “t1e7e" of3 dofdw yd ebordtem eetdy sta ovedT sew! 
af? gaidoeteb sviovat saed3 to owT ,Inomrsent atl? x — 
yaseb eft mott gotiivass eslotiisq begrsdo bas eyat~ammqeg: baast 
gaieqtise mt bemzo2 .vlavidssques bavodnay baa bawod , ‘hei 
moo bes aevplofoes anteluq masd ayolqas bodiem brtds sat ts . 
.ehorsuen sift griouborq notiseex to sqy3 si2 to snsbase 
z0% teazved sid otno beeluq at gelotizeq Yo mead potbedaod 0 
Yo 398% & 35 ,bnoosacnsc ono ylieotays ,aat? Yo abot 194 _ 
[fgteb ni bodlioesh at mesd odt seluq 02 boaw bods aft Fe 
saveood vitetrd sre cxnesnean {Iiw 3ud asoneTster avo r ‘; 


Yo asgnometassom vissimilera omoe mt beyolqus eaw sup hiv 


V aon: ae vig 
booubat fangte s mort bevirsb sts esalug somer9} 1) 
ob eds gntqoewe ok bsev sgsilov Isbloeuats 3 13 yd goo. 


evtie "seqanor" oft seoms baad tes: 2 aD ai : meV 


S| oat 60 ie en 


Le a ai aor a 
ay sy © a eee . 4% we 2 ae) ae 


-~ 27 = 


(see Figure 2-1), This signal is shaped and squared to produce one 
positive, 150 nanosecond wide, "start" pulse for each rf cycle, 
There are two current bursts produced by the pulsing system per rf 
cycle, only one of which can be used (providing two timing pulses 
per rf cycle would be unsatisfactory because the two beam pulses 
form separate spots on the target, have different flight times to 
the target, resulting in increased time dispersion, and probably 
contain different amounts of current). The unwanted pulse can be 
eliminated electronically by the Alternate Pulse Eliminator, or 
mechanically by means of an adjustable slit system in the drift 
tube immediately following the analyzing magnet (utilizing the fact 
that the two pulses are physically separated in space). A 300 
nanosecond delay line ireleot AS=48) is used to obtain the correct 
phasing between the "start" pulse and the beam pulse. 

Another method of obtaining the "start" pulse using the pulsed 
beam is to sense the beam pulses at the target, or at the Faraday 
cup, and convert them into the zero-time reference by means of the 
electronics described by Alexander (Al 63). This method has the 
particular advantage of eliminating the variation in the flight 
time of the beam burst in its travel down the drift tube from the 
point of formation approximately 6 meters distant from the target, 
This dispersion results from the variation in the incident beam 
energy, which is approximately +3 keV. 

In the method of obtaining the "timing" pulses by detecting 


the associated charged particles from the decay of unbound states 
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formed in stripping reactions, a solid state detector, located close 
to the target, is employed, advantage being taken of the precise 
timing (Al 63) that one can achieve with this type of detector, This 
is the method used in the present studies of the alpha particle decay 
of gee" formed by the st (a nyc!** reaction, The fast electronic 
circuitry required for fast timing with solid state counters is 
described in detail by Alexander (Al 63) and will be mentioned only 
briefly in the following discussion of the electronics involved in 
the neutron spectrometer. 

Figure 2-5 shows a block diagram of the complete neutron spec~ 
trometer as used in the. present experiments, as well as the charged 
particle spectrometer employing two solid state detectors. The time 
converter "start" pulses are initiated by the particles detected by 
the solid state counter labelled Ops which is a Nuclear Diodes” low 
resistance, surface barrier device, type PL=2=50=10, operated at a 
bias of -180 volts, The charge sensitive preamplifier, which employs 
two 8056 nuvistor triodes in the input stage, integrates the detector 
pulses, preparing them for double delay line clipping in the following 
amplifier, The bi=-phased output of the amplifier (each phase having 
a width of 100 nanoseconds) is further amplified by a Hewlett Packard 
wide band distributed amplifier, model 460AR, before being fed into 
the cross-over-pick-off circuit. This circuit generates the reference 


pulses for the time converter, and their shape is shown in the figure. 


The front edge of the positive portion of the pulse is timed from the 
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"point" (in time) when the input pulse crosses the zero axis, Its 
width is determined by a 150 nanosecond shorted delay line. The 
negative going portion of the pulse is retained to keep to a minimum 
base line shifts which would result from high counting rates, With 
the exception of the nuvistor tubes used in the preamplifier, and 
the distributed amplifier, transistors were used throughout in the 
foregoing circuitry. 

The neutrons were detected by NE 213 liquid scintillator 
contained in a quartz cell (6 in. x 3 in. x 1 in.) which is viewed 
by two photomultiplier tubes. The XP 1040 feeds a limiting amplifier 
with delay line clipping that produces the 150 nanosecond "stop" 
pulses for the 6BN6 type time-to-amplitude converter, It also drives 
a linear preamplifier which feeds the neutron side channel formed by 
a fast amplifier and a single channel pulse height selector. Thus, crude 
pulse height selection on the neutron spectrum is provided. This side 
channel is usually set to eliminate those pulses which are too small 
to produce properly limited pulses for the time converter, 

The RCA 7046 photomultiplier feeds a pulse shape discriminating 
circuit which distinguishes events in the scintillator initiated by 
neutrons (recoil protons) from those initiated by gamma-rays. The 
output from this circuit is illustrated in Figure 2-5. The gamma-ray 
eliminating side channel, consisting of an amplifier (two x5 general 
purpose, transistorized amplifiers each having a delay of 0.25 


microseconds) and a single channel pulse height selector, can be set 
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to almost completely eliminate random counts resulting from gamma-ray 
detection in the liquid scintillator. 

In addition to requiring the simultaneous occurrence of the 
positive pulses from the two above mentioned side channels, the 
kicksorter gate generator requires the absence of pulses from the 
negative-time-eliminator, This circuit produces a negative pulse 
whenever, by chance, the "stop" pulse precedes the "start" pulse 


into the time-to-amplitude converter. 


The factors which influence the timing resolution of the 


spectrometer are listed and discussed by Neilson et al (loc, cit.) 


and Riley (Ri 62), and need not be repeated here. Suffice it to 
say that for alpha particles detected in coincidence with gamma 
radiation, Alexander (Al 63) found the timing resolution to be 2 
nanoseconds, He used the tae CT hea reaction to produce 2 MeV 
alpha particles in coincidence with 6.13 MeV gamma radiation, 

The energy resolution is related to the timing resolution by 
the expression AE/E = 2At/t, where E and t are the energy and the 
flight time of the neutron respectively, and At is the time resolu~ 
tion (At = 2 nanoseconds). In the present experiment, the maximum 
energy resolution could not be realized because the low counting 


rates limited the flight path to one meter, 


VI. The Charged Particle Spectrometer 


As mentioned in the introduction to this chapter, two types of 
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experiment were performed in this investigation of the alpha decay of 


12 ; 
C "3; one using a single solid state counter and the other using two 


solid state counters, The complete electronic arrangement for the 
latter type is shown in Figure 2-5, The detector Ot, is used to detect 
one or two of the alpha particles from the decay of a given state in 


2 
c} while Op is used to detect the other remaining particle. Os is 


another Nuclear Diodes surface barrier detector, similar to Op and 
operated at the same bias, It feeds a nuvistor preamplifier, iden- 
tical to the one already mentioned, which in turn drives a 1 micro- 
second double delay line linear amplifier (Go 60a). The signals from 
this amplifier are gated through a linear gating unit (Go 60b) by 


pulses obtained from a,, which were amplified in the same way as those 


yi 


from Ops This unit also uses the zero cross-over technique to 
determine the delay between the two input pulses A and B, Thus, a 
resolving time of 35 nanoseconds is realized with this instrument. 
The linear gate output, which then corresponds to charged particles 
in coincidence with those detected by the "timing" detector, is fed 
to the horizontal input of a TMC, 1024 channel, two parameter kick- 
sorter operating in the 32 x 32 mode. The vertical input of the 
pulse height analyzer takes its input from the time-to-amplitude 
converter so that the coincident particle spectra of the horizontal 
input are sorted out according to neutron groups to which they belong 
ké 


or, in other words, according to the level of excitation of C 


from which they originate. 


The Oy pulses are self-gated through a second linear gate unit 
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and fed to a 256 channel TMC kicksorter so that the operation of this 
detector could be monitored continuously. 
When performing the experiments requiring only one solid state 


detector, a, and its associated electronics were removed. Pulses 


E 
from the "timing" detector were self-gated through the linear gate 
unit and fed to the horizontal input of the two-parameter analyzer. 
The neutron time spectrum again provided the vertical input. The 
kicksorter was, of course, as in the former case, gated by the side 
channels already discussed, 

The 256 channel pulse-height analyzer was used to obtain neutron-= 
gated particle spectra, and neutron time-of-flight spectra to compare 


with that obtained when using the beam pulsing system and isotopically 


yes * 
separated B targets . 
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* Obtained from A.E.R.E., Harwell, England 
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CHAPTER III 
EXPERIMENTAL MEASUREMENTS 


I. Introduction 


This chapter deals with the measurements made in the investigation 
of the alpha-particle decay of the 12,71 and the 11,83 MeV states in 
rae formed by the aie | ht anal reaction, 

Section II describes preliminary measurements of the neutron 
time-of-flight spectra obtained using the pulsed beam technique, and 
an isotopically separated ta target as well as the self-supporting, 
natural boron target used in the subsequent experiments. 

The measurements made using a single solid state detector, from 
which were obtained the "timing" pulses and the charged particle 
spectra, are described in section III. One and two parameter data 
are presented and discussed, 

Section IV is concerned with the experiments in which two solid 
state counters are employed in an effort to detect simultaneously two 
or all three of the alpha particles from the decay of a given level 
of a Two parameter data which accentuate the 12.71 and the 11,83 
MeV levels in c)? are presented. 


The conclusions that can be drawn from the experimental data are 


summarized in the final section, V. 
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Figure 3-4, Plan view of detector geometries. 
and oc are '!Nuclear Diodes! surface barrier 
detectors. 


ll | 


eo eer Aa ages Sect a re oe eT eee 


eS TO Se see eS ae 


aot teat) 
rg Ioet eo 


“houses 
so2%5o78D 


a ae sree ae eee 
fe .2@3i739mos9. 107S5I8b Fore | 
salased sostiwe ‘esbold wesrel 


a ah 


II. Pulsed Beam Neutron Time-of-Flight Spectra 


A time-of-flight spectrum for the neutrons produced in the 
e (d,nc-- reaction was obtained using the pulsed beam technique 
and an isotopically pure pit target , Lest nie fen thick, to provide 
a "standard" with which to compare the neutron spectra obtained 
in the experiments using the self-supporting, natural boron targets, 
and associated particle timing. Such a spectrum is shown in Figure 
3-1, It was taken at a deuteron bombarding energy of 1 MeV with the 
neutron detector at an angle of 25° with respect to the beam, and 
2.8 meters distant from the target. The gamma-ray-eliminating side 
channel and the neutron side channel were set to cut off pulses 
produced in the liquid scintillator by a RdTh gamma-ray source, 
In the case of the neutron side channel, this is equivalent to 
setting the neutron threshold at about 600 keV. All the levels 
in oe up to 12,71 MeV excitation, with the exception of the some- 
times postulated 10,1 MeV state, are evident in this spectrum, It 
is also apparent that the levels above the alpha-particle breakup 
threshold, 7,28 MeV, are riding on what could be a neutron continuum 


La 
produced, perhaps, by the simultaneous three-body breakup of C 


The three-body breakup (a neutron and two alpha particles) of Be” 
produced by the B*(a.0)Be reaction may also contribute to this 
continuum, On the other hand, the hump in the time spectrum may not 
be due to a continuum at all, but may be the result of the overlapping 
eee nnn 


* Separated tag target obtained from A.E.R.E., Harwell, England 
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levels and a very broad level at 10.1 MeV (Ni 62). It is one of the 
aims of this research to help establish the origin of the hump in the 
time spectrum, 

Figure 3-2 shows the neutron time-of-flight spectrum taken using 
the self-supporting, natural boron target, Running conditions were 
the same as those mentioned above for the separated target except for 
a slight difference in the gain of the kicksorter, The dashed curve 


near channel 100 indicates the contribution of the pl? in this target 


to the time spectrum, It was obtained by using a self-supporting ae 
target, and by applying the appropriate normalizing factors. The 
effect of the a on the time spectrum from the natural boron target 
is completely eliminated when the associated charged particles are 
used to generate the zero-time reference for the time-to-amplitude 
converter, unless the solid state detector is able to "see" the ae 
recoil nucleus. This possibility was removed by using a 10 yin. 
nickel foil in front of the detector, 

The calibration of the generating voltmeter of the Van de 
Graaff accelerator was obtained from the measurements made by T.K. 
Alexander of this laboratory on the thickness of the boron target as 
described in Chapter II, section IV, The thick target gamma-ray 
yield from the 873 keV resonance in the reed eT ees tad reaction was 
measured as a function of the generating voltmeter reading, Curve A 
of Figure 3-3 shows the yield curve obtained. The resonance energy 
of 873 keV is seen to correspond to a meter reading of 858 keV. This 


discrepancy was taken into account when setting the bombarding energy 


of further experiments. 
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(see text). 
(Courtesy of T.K. Alexander) 


macau! 


SR SS Te Oe Sa RES SE RS I SS SS IEEE EE ee ee 


lat 


mearese es so eee ems 


O02 use 
sibeen a eow Fe 
pal 2 25viU5 bf iy ty 
-*ig2 Bae Asixsiditss ese 
InemsiWesam 229715 As septs) nope 


(FebnsxslA WaT 
> 


- 36 =- 


TII. Measurements Using One Solid State Detector 


(a) Experimental Arrangement and Running Conditions 


The detector geometry used in the experiments employing a 
single solid state detector to provide both the timing pulses and 
the charged particle spectra is shown in Figure 3-4a, The angles 
of the detectors and the target are measured from the direction of 
the incident beam, and are positive as indicated by the arrows. 
Because these experiments did not require the best possible reso- 
lution from the solid state detector, which was less than one 
percent with good geometry and an Ortes low noise amplifying 
system, a large solid angle was used in an effort to achieve a 
reasonable (n,a) coincidence rate with a neutron flight path of 
one meter. The detector, Ap» was shielded by a brass collimator 
with a 5/16 in. diameter hole in the center, and was positioned 
31/32 in. from the beam spot on the target. This corresponds to 
a solid angle subtended by the detector at the target of 8.2 x 
id steradians, or a half angle of 9.1 degrees, 

The neutron flight path for most of the runs was 1,15 meters 
although one time spectrum was taken at 1.42 meters. The counting 
rate at this distance however was almost intolerably low, The 
energy resolution at 1,15 meters and a bombarding energy of 1.5 
MeV for the lowest energy neutron group of interest (that is, those 
leading to the 12.71 MeV state in chi was 8.5 percent. 


cqmnepeneeesrepeernttpenemnenetprtpeentasnet cae EE SE SLI 


* Oak Ridge Technical Enterprises, Oak Ridge, Tenn. Model 101-201 
low noise preamplifier-amplifier system. 
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Since the 12.71 and the 11.83 MeV states in a are of particular 
interest in this work, the neutron detector was set at angles corres- 
ponding to the stripping peaks of these levels. For the 12.71 MeV 
level, the dueteron stripping peak occurs at approximately 45 degrees 
in the laboratory system (the angular momentum of the captured proton 
being almost certainly _ = 1), and for the 11,83 MeV level it occurs 
at an angle of zero degrees, P being zero for this case (Sa 62). 

A series of runs was carried out with the neutron detector at each 
of these angles. 

A deuteron bombarding energy of 1.5 MeV was used throughout 
in the following experiments, and the beam current was kept below 
0.1 wamp. at which the total counting rate in the solid state detec- 
tor was of the order of 20,000 per second. 

The beam was collected from the target holder and a Faraday 
cup, both biased at +300 volts to suppress secondary electron emission, 


* 
and measured by an Elcor current integrator, model A309A. 
(b) Setting-Up Experiments 


In order to check that the system was operating properly, some 
preliminary, particle-gated neutron time-of-flight spectra were 
obtained. Figure 3-5 shows the spectrum recorded by a 256 channel, 
TMC pulse-height analyzer for y = 90°, oT = 45°, and oF =z 45°: 

Sn = 1,15 meters, The peaks were identified by comparing their 
separation in time with that calculated from the energies of the 
ee 


* Elcor Inc., Falls Church, Virginia 
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neutrons leading to the indicated states in e*. The effectsof the 


alpha-particle flight times to the solid state detector on the time 
separation of the neutron peaks were negligible. The kinematics 
were calculated using the IBM 1620 computer program Nuclear Reac- 
tions I, fully described in Appendix I. The four states in po 
above the alpha=-particle breakup threshold which appeared in the 
pulsed beam spectra, Figures 3-1 and 3-2, are also apparent in this 
spectrum, In addition, however, there is a peak which corresponds 
to the 4.43 MeV state of gh which is well below the breakup 
threshold. That this peak was in fact due to the neutrons leading 
to the 4.43 MeV state, having been "timed" by the recoiling 
B"(4.03) nuclei was shown in the following way, Time-of-flight 
spectra were recorded with the solid state detector (1) in line 
with these recoiling nuclei (y = 105°) and (2) at an angle (y = 80°) 
such that they were not within the solid angle of the detector. 
Figure 3-6 shows the results of this experiment. Curve (1) was 
obtained with 35 uc of charge and shows a very prominent peak at 
the 4,43 position, Curve (2) on the other hand required 110 uc of 
charge and there is no evidence at all of a peak at the 4.43 
position, That the co recoil nuclei have little difficulty in 
getting out of the target was further substantiated by T.K. Alexander 
of this laboratory who detected them in coincidence with 4.43 MeV 
gamma radiation (Al 63a). 

To eliminate this added complication to the neutron time 


spectrum, a 10 pin. nickel foil was placed in front of the solid 
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state detector, thus preventing these recoiling heavy nuclei from 
starting the time converter. Unless otherwise stated, all sub- 
sequent runs were made with this foil in place, 

Figure 3-7 shows the time-of-flight spectrum taken at a 
neutron flight path of 1,42 meters; y = 90°, oe = 45°, and 
Ps = 45°, The resolution here is not significantly better than 
that obtained for i = 1,15 meters and the running time was over 
three times as long (Figure 3-5) to obtain the same statistics, 
The running time for the spectrum of Figure 3-7 was 116 minutes 
which indicates an almost intolerably low coincidence counting 


rate, In view of the above, the flight path used in the remaining 


measurements was set at 1.15 meters, 
(c) One=-Parameter Measurements at 9. = 45° 


Measurements with the neutron detector at 45° were made 
in order to accentuate the 12,71 MeV state in an as mentioned 
in part (a) of this section, Two 256 channel, TMC pulse-height 
analyzers were used to simultaneously record the charged-particle- 
gated neutron time-of-flight spectra and the neutron-gated charged~ 
particle spectra as functions of the solid state detector angle, 
y. Both kicksorters were gated by the coincidence circuit gate 
generator (shown in Figure 2-5) and the side channels were set 
as previously described. 

Spectra were taken for thirteen positions (Yq) of the solid 


state detector from 60° to 265°. Direct charged-particle spectra 
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for each setting were also recorded, Since the variations in the 
character of these curves from setting to setting were not rapid, 
only the data for four positions representative of the trends are 
presented, These data are shown in Figures 3-8 to 3-11 for 

y = 90°, 130°, 240°, and 265°, The peaks in the direct spectra 
were identified with the aid of kinematical calculations made for 
all deuteron reactions that could occur in the target material and 
its possible contaminants. The calculations were performed by the 
University of Alberta IBM 1620 computer using the program, 

Nuclear Reactions I, described in Appendix I, The ad alpha- 
particle source in the target chamber was used to provide the 
energy calibration of the kicksorter. 

Computer calculations have also been performed to determine 
the energies of the alpha groups which would result from the 
decay of the states in e through the ground and first excited 
states of he” These calculations are fully described in Appendix 
I, under the heading Nuclear Reactions II, and the results for the 
decay of the 12.71, 11.83, and the 9.64 MeV levels of C’* are 
given there, 

If the 12,71 MeV state decays to the ground state of be’, a 
rather sharp group of alpha particles is expected in the neutron~ 
gated charged-particle spectra, These spectra, of which Figures 
3-8a to 3-lla are typical, show no evidence of this group. This 
is understandable on the basis of the tentative spin-parity 


assignment of \* (Aj 62) for the 12,71 MeV state, making alpha 
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decay to the of ground state of Be® forbidden, Hence, these data 
provide very strong support for the iF assignment for this level 
of okt; 

Also, these spectra do not show peaks (where indicated in the 
figures) which would correspond to the alpha group from the 11,83 
MeV level of otf leading to the Be® ground state, Thus, if the 
11.83 state does decay sequentially through Be®) it would seem 
that it does so predominantly via the 2.9 MeV state. There is no 
selection rule forbidding the unobserved transition since there 
is little doubt about the 1 spin-parity assignment for this level 
of ces (Sa 62). 

The two peaks which occur near the end of the charged- 
particle coincidence spectra, as for example in Figure 3-9a at 
3.4 and 4.35 MeV, are attributed to alpha particles from the 
pt cataype? (3,04 and 2.43 MeV) reactions respectively in 
coincidence with neutrons from the subsequent breakup of the Be? 
states, The 3.4 MeV peak may also have contributions from random 
coincidences between neutrons and protons from the eo? (aot? tend) 
reaction, 

With reference to Figure 3-8a, the peak at 1 MeV and the 
hump at 1.25 MeV are due to random coincidences between neutrons 
and scattered deuterons and the very intense proton group from 
the Bt a p)Bt* (0.95) reaction respectively. The peak near 2 
MeV consists, in part, of random events involving protons from 


11 16 17 
the Bt (a,p)p!2 (gna), BLO (d,p)B (7.99), and 0 (d,p)O° * (0.87) 
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reactions. It also contains alpha particles from the decay of the 
9.64 MeV state of gi? to the ground state of pes a well established 
transition (see Chapter I), and those from the decay of the 12.71 
MeV state of ole to the 2.9 MeV state of Be”) assuming that this 
transition takes place, The valley at 1.5 MeV fills in, and the 
peak at 2 MeV decreases in height as y increases throughout its 
range; this is in keeping with the kinematic shift in energy of 

the alpha-particle groups from the 9.64 and 12.71 MeV levels of 
ae", The overlapping of these two alpha-particle groups points out 
the need for a two parameter analysis of this case to sort out the 
charged=particle spectra according to the neutron groups with which 
they are associated or, in other words, according to the levels of 
gh? from which they originate. Such measurements were carried out 
and will be reported in part (d) of this section. 

Figures 3-8b, 9b, 10b, and 11b show the charged-particle-gated 
neutron time-of-flight spectra corresponding to the charged-particle 
spectra discussed above. One of the features of these curves is 
the variation in the intensity of the neutrons leading to the 12.71 
MeV state of os with the position of the solid state detector. It 
is clear that there are more alpha particles originating from this 
state on the side of the deuteron beam opposite the neutron detector 
in the general direction of the recoiling o'? nuclei, Now if this 
level decays by a direct three-body process, there would be a high 
probability of equal sharing of the disintegration energy between 


the three particles as indicated from the three-body distributions 
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shown in Figure 3-15, The total disintegration energy being 5.43 
MeV, this means that each particle would have about 1.81 MeV of 
energy. Therefore, the recoiling a nuclei, which have an energy 
of only 0,134 MeV (at an angle of 60° from the deuteron beam), 
would be expected to have little or no effect on the number of 
alpha particles emitted in a given direction in the reaction plane, 
and hence the angular distribution should be almost isotropic. 

On the other hand, if the decay of this state proceeds through 
he" (in particular, the 2.9 MeV level of Be°) , it can be argued 
that one should expect the angular distribution of the alpha 
particles to be peaked about the center of mass direction of the 
recoiling e nuclei, Assuming isotropic breakup of the fee state 
formed (2.9 MeV), this angular distribution would be reflected in 
the neutron time-of-flight spectra, 

The angular distribution of the alpha particles from the 
decay of the 12.71 MeV state to the 2.9 MeV state of Bae can be 
calculated by assuming that the carbon level is an intermediate 
stage of the reaction, a in gipe ook This assumption is 
reasonably well justified because the stripping pattern for the 
neutrons leading to the 12.71 MeV state is predominantly Butler 
with the angular momentum of the captured proton equal to one 
(Zp = 1) (Sa 62). The angular distribution will then be symmet- 
rical about the direction of the captured proton; that is, about 
the direction of the recoiling e nucleus (To 61). Taking the 


ground state of pit as having a J” of a/2 . there are two possible 
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values of the incoming-channel spin; Sy = 1 and Sy = 2, The 


incoming orbital momentum is a = 1 (Sa 62), the intermediate cl? 
+ 
state is taken as al level, and the out-going channel spin and 
the out-going orbital momentum are both 2, based on the a" assign= 
8 
ment of the Be (2.9) state, From the tables of Sharp et al (Sh 54), 
and using the above angular momentum quantum numbers, the angular 


distribution functions were found to be, for S, = 1, 


1 
W,(¥) ~ 14 “ P, (Cos) , 


and for Ss) = 2, 


1 
W, Cv) v~ 1-45 P, (Cosy), 


where » is the angle of emission of the alpha particles in the 
coordinate system in which the a is at rest, P, (Cos) is the 
usual Legendre polynomial of the second degree. These functions 
are shown in Figure 3-l2a plotted against ~ and the corresponding 
laboratory angle y. Also shown in the figure are experimental 
points representing the peak heights of the co (12.71) neutron 
group in the thirteen measured particle-gated time-of-flight 
spectra, These points have been normalized to 500 uc of charge, 
to a target angle of 45°, and to 1.5 at y ® 90°; they have also 
been corrected for background and for the solid angle transforma- 
tion from the laboratory system to the center of mass system of 


the on nucleus (see Appendix I, Nuclear Reactions I, where the 
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kinematics have been described in some detail). the shape of the 


experimental distribution compares reasonably well with W The 


1° 
shift in the symmetry axis, approximately 40° in the center of mass 
system, from that defined by the recoiling gts nuclei could be a 
manifestation of the distortion effects taking place in the forma- 
tion of the carbon state (To 61). That such effects are present in 
this case is indicated by the rise in the neutron angular distribu- 
tion at back angles (Sa 62). 

On the basis of the agreement of experimental distribution of 
Figure 3-12a with Wy and not with an isotropic distribution 
(which would be a straight, horizontal line in the figure), it 
would appear more likely than not that the 12.71 MeV state of ge 
decays sequentially through the 2.9 MeV level of Se’, It would 
also seem that Sy = 1 is the dominant incoming channel spin for 
the reaction, 

Because the peaks in the time-of-flight spectra of the group 
corresponding to the 11.83 MeV state in qt are not very prominent, 
an angular distribution based on the peak heights would hardly be 
justified, Nevertheless, the experimental points, normalized and 
corrected for solid angle, are presented in Figure 3-12b in order 
to get a crude idea of the distribution of alpha particles origi- 
nating from the 11,83 MeV state. The predicted distributions for 
the alpha particles leading to both the ground and first excited 


states of Bee are isotropic, assuming that pure Butler stripping 


is the dominant process in the formation of the carbon state, with 
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‘. = 0 (Sa 62), and that the He" is formed by a a (ocnlbes reaction, 

Two of the three three-body energy distribution functions for 
this case, shown in Figure 3-18 (curves a and b) and discussed in 
detail in part (f) of this section, predict equal sharing of the dis- 
integration energy. On this basis, a three-body breakup of this 
state would also be isotropic since the energy of the recoiling ch 
nuclei (.153 MeV) would not significantly alter the distribution, 

The very crude experimental distribution of Figure 3-12b could 
be interpreted as being isotropic. This only confirms what would be 
expected irrespective of the mode of decay. 

The experimental angular distribution of the alpha particles 
from the 9.64 MeV state of ge to the ground state of bee is shown 
in Figure 3-13, These points (open circles) correspond to the 
peak heights of the ct (9,64) neutron group in the neutron time= 
of-flight spectra, corrected as before for background and the solid 
angle transformation, They are normalized to 1,58 at y = 90°, and 
are connected by the dot-dashed curve for visual aid. The three 
points enclosed in triangles were obtained from two-parameter data 
and will be discussed in part (d) of this section, The solid curves 
shown in Figure 3-13 are the distributions predicted assuming a 
proton capture reaction in the direction of the (9,64) recoil, 
and were determined using the following angular momentum and parity 


a 12 a 
assignments (Aj 62): * cv ae A. = 2 (Sa 62); C (9.64), 3 ; 


pe°(0), bh. W. is the distribution for the incoming channel spin, 
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Sis equal to one, and is given by 
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W,) ~ 1+ 8/7 P. (Cosy) + 6/7 P, (Cosy) . 


The distribution for the other possible incoming channel spin, 


Sy = 2, is 


W, (v) ~ 1+ 2/7 P, (Cosy) - 9/7 P, (Cosy). 


We, is W, + »25W) « Ws and W, are the distributions one gets for 
the two values of the incoming-channel spins using a l- assign- 
ment for the ct (9,64) state (Aj 59), with everything else the 


same as above. Explicitly, 
, 
W,@) ~ 1+ P,, (Cos) 
and 
' = 
Wy) v1 P,, (Cosy) ° 


It is evident that no combination of Wy and W, could be 
found that would reproduce the general trend of the experimental 
curve, Hence, the 1? assignment is ruled out. On the other hand, 
it is equally evident that Wo plus a small amount of Wy rather 
closely resembles the experimental distribution, A smearing or 
flattening of the distribution and a shift in the symmetry axis 
are explicable considering the inability of the Born-approximation 


stripping theory alone to account for the increase in the neutron 


intensity at back angles (Sa 63). Therefore, the experimental 
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distribution of Figure 3-13 is considered to be strong evidence in 


favor ot the 3. assignment for the 9,64 MeV state of 6, 


(d) Two Parameter Measurements at 8 = 45° 
n 


Two parameter analysis of the neutron time-of-flight spectra 
and the neutron-gated charged-particle spectra were carried out to 
obtain energy distributions of the alpha particles originating from 
the levels of cl? fed by the (d,n) stripping reaction and, where 
possible, to obtain angular distributions of these particles, 

With the neutron detector fixed at 45° in the laboratory system, 

the 12,71 MeV state is again the accentuated level of a” as 
discussed earlier, The experimental arrangement was as described 

in part (a) of this section with the neutron flight path set at 

1,15 meters, The 1024 channel, two parameter, TMC pulse-height 
analyzer operating in the 32 x 32 mode was used. Data were recorded 
for three positions of the solid state detector: y = 70°, 90°, and 
120°. For angles greater than 120°, the (n,a) coincidence rate was 
too low to be practicable. As it was, the running time for each of 
the three mentioned angles was approximately three hours, 

An isometric view of the two-parameter data obtained for 
y = 90° is shown in Figure 3-14, The charged-particle-gated neutron 
time-of-flight spectrum is displayed on the y-axis, and the neutron- 
gated charged-particle spectra are displayed on the x~axes. The 


charged-particle spectrum at y = 0 corresponds to all the charged 


particles which are in coincidence with all the neutrons of the time 
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spectrum (as in Figure 3-8a), and the others correspond to only those 
charged particles which are in coincidence with neutrons in the 
adjacent part of the neutron time-of-flight spectrum. 

Of particular interest here is the spectrum of the charged 
particles which are in coincidence with the neutrons leading to the 
12,71 MeV state of Gr, This spectrum is primarily contained in the 
one labelled no, 3 in the figure and, to a lesser extent, in spectrum 
no. 4, 

Figure 3-15 shows spectrum no. 3 redrawn in greater detail, The 
large peak at the low energy end of the spectrum is the result of 
random coincidences between the neutrons and scattered deuterons. 

The two broad peaks at 2 and 3,1 MeV are only consistent with the 
interpretation that they are due to alpha particles from the sequen- 
tial decay of the 12,71 MeV state through the 2.9 MeV state of Be’, 
with the 3 MeV peak "i ~ ,8 MeV) resulting from the breakup of the 
Bae state, and the 2 MeV peak Soe Ay ~ .7 MeV) from the alphas 

leading to the Be® state, 

Also shown in Figure 3-15 are three curves a, b, and c which 
represent the type of energy distributions one might expect if the 
12.71 MeV level disintegrates by a direct three-body process, The 
three curves are neither normalized to each other nor to the 
experimental spectrum as they are presented only to show the general 
shapes of the three-body distributions. 


Curve 'a' is the classical distribution (Ne 58) given bv 
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where € = E/E ax’ E is the center of mass energy of any one of the 
particles and ee ee is the maximum kinetic energy that it can have 
which, in this case, is 2/3 Q (Q, = the disintegration energy = 
5.43 MeV. The kinetic energy of the c1* (12,71) nucleus has been 
neglected as it is not expected to significantly alter the shape of 
the distribution. 

Curve 'b' is the classical distribution multiplied by the 
product of the three penetrability factors for the alpha particles. 
It was felt that such a function would at least illustrate the 
general nature of the distribution with centrifugal and coulomb 
barriers taken into account, Explicitly curve 'b' is given by 


C1) <2) G3) 
f,(€1) = f fe Py P P 


where P denotes the penetrability of the abn particle emitted 


(n) 
Q 
with orbital angular momentum 2, The & values assigned were based 
on a ie spin-parity assignment for the 12.71 MeV state, Two alphas 
with one unit of orbital momentum or one with two units (the others 
having zero units) are the only possibilities for this case, and 
of these two, the former would be more probable than the latter. 
The assumption was made that for a given value of ay the other 
two alpha particles shared equally the remaining energy. 

Curve 'c' shows the distribution one would expect using the 


function given by Delves (De 60) which takes into account the 


angular momentum barrier and the three-body impact parameter. In 
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general, this function is 


hy + 1/2 he + 1/2 
f(e) =e (1 - €) ; 


where hy + 7) = J, the angular momentum of the decaying state, In 


particular, for this case the most probable value of 2, and 2. is 


i 2 


again l, so that 
f (e) = 3/2 - gy rte ‘ 


All three distributions indicate that there is a high probability 
that there would be an equal sharing of the disintegration energy 
among the three particles (and, hence, that the angular distribu- 
tion of the particles would be isotropic). Such a sharing is 
clearly inconsistent with the experimental energy distribution, 

The three points enclosed by triangles in Figure 3-12a were 
obtained from the peak heights of the 2 MeV group of Figure 15 and 
from the corresponding groups in the data for the other two angles, 
y = 70° and 120°, They were normalized to 1,5 at 90° as well as 
corrected for the solid angle factor. They are in very good agree- 
ment with the data taken from the single parameter, time-of-flight 
spectra, 

Very little can be said about the 11,83 MeV state of ha 
from this two parameter data. If this level does decay through 


Be°(2.9), one would expect broad peaks in the regions of 1,5 MeV 


and 2-2,5 MeV, There are traces of such peaks in the charged- 
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particle spectrum (Figure 3-14, no. 5) associated with the 11,83 
MeV level but the statistics are very bad since the neutron detector 
is well off the stripping peak for this level, 

The decay of the 9.64 MeV eb state to the ground state of ie 
shows up very well in spectrum no. 8 of Figure 3-14, The effect of 
this level is also evident in spectra no, 7 and 9. Similar peaks 
were obtained in the data for the other two positions of the solid 
state detector (70° and 120°), Measurement of the area under the 
peaks due to the 9,64 MeV level for the three detector positions 
provided the three points of Figure 3-13 enclosed in triangles, 
These agree very well with the points of the angular distribution 
obtained from the single parameter time-of-flight spectra discussed 


in part (c)i. 
(e) One-Parameter Measurements at a = 0° 


With the apparatus set up as described in part (c) of this 
section, one-parameter neutron time-of-flight spectra and neutron- 
gated charged-particle spectra were obtained with the neutron 
detector on the stripping peak of the 11.83 MeV level of er zero 
degrees, for five settings of the solid state detector: y = 70°, 
90°, 100°, 120°, and 140°, The 140° data are presented in Figure 
3-16, 

The angular distribution of the alpha particles from the 
decay of the 11.83 MeV state, as implied by the normalized and 


solid-angle-corrected neutron peaks corresponding to this state, 
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is isotropic within the statistics of the points. The corrected 


data are as follows: 


Corrected neutron 


y Counts 
70° 167 + 13 
90° 165. @) 13 

100° 179 # 14 
120° 173 + 13 
140° 194 + 14 


There are four peaks in the neutron-gated charged-particle 
spectrum of Figure 3-l6a. The two at 3,35 and 2.75 MeV are due to 
coincidences between the primary alpha particles and the neutrons 
produced in the reactions: ora aye” (3.04 and 2.43 respectively), 
Re Ss n+ ie”, The protons from the 6 taene Gand) reaction would 
also contribute some random coincidences to the 3.04 MeV peak, The 
peak at 1.5 MeV is due to alpha particles from the 12,71 MeV state 
of ate to the 2.9 MeV state of i. and that at approximately 2 Mev 
can only be accounted for by alpha particles from the breakup of the 
2.9 MeV state of Be’, fed by any or all of the levels in he at 
10,84, 11.83, and 12.71 MeV. Alpha particles from the 10.84 MeV 
state to the Ha’ ground state, if this transition occurs, would 
contribute to the hump at 2,5 MeV, If the 11.83 MeV level decays 
through Be°(2,9), a broad peak at 1 MeV would be present; such a 
group would, however, be lost in the large peak resulting from 


random events. There is no evidence in any of the five charged- 
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particle coincidence spectra for peaks corresponding to transitions 
from the 11,83 or the 12.71 MeV states to the Be” ground state, In 


Figure 3-l6a, these would occur at 2,9 and 3.5 MeV, respectively. 
(f) Two-Parameter Measurements at eT = 0° 


On the basis of the single-parameter spectra it can only be 
said that the 12,71 MeV state of Er does not decay to the ground 
state of ta”, a conclusion which is in agreement with the measure- 
ments made at a = 45°, Nothing, however, can be said regarding 
the disintegration of the 11.83 MeV level because of the complexity 
of these spectra, It was hoped that a two-parameter analysis of the 
neutron and charged-particle coincidence spectra would provide a 
definite indication of the decay mode of this state, and perhaps of 
the 10.84 MeV level as well since its stripping peak is also at 0° 
(Sa 62). The experimental arrangement was as described in part (d) 
of this section, with the neutron detector at O° and the solid 
state detector at 90° and 140°, 

An isometric drawing of the data obtained for y = 140° is 
shown in Figure 3-17. The data taken at y = 90° has the same 
general features so that the following remarks apply to both cases. 

Spectrum no. 5 shows the charged particles in coincidence 
with the neutrons leading to the 11.83 MeV state; it is redrawn 
in Figure 3-18 to a larger scale to facilitate closer examination. 
The energy scale was obtained by comparing the "single-parameter' 
spectrum, no. O of Figure 3-17, with that of Figure 3-16a and with 


a direct spectrum taken at 140°, The various alpha particles 
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resulting from the sequential decay of the 11,83 MeV state through 
Be® are shown in the energy level diagram in Figure 3-18, Their 
positions on the energy axis of the spectrum are also indicated, 
allowance having been made for the energy losses in the 10 yin. 
nickel foil in front of the solid state detector. Peaks corres- 
ponding to transitions to the ground state (246) and to the breakup 
of the 2,9 MeV state of hee are evident. The latter peak, at 
approximately 1.5 MeV, no doubt contains alpha particles from the 
2.9 MeV level of ie” fed from the 12.71 MeV level of ras (as this 
level overlaps the 11.83 MeV level in the neutron time-of-flight 
spectrum), but must be due primarily to the alpha particles, 


a from the 2.9 MeV level fed by the 11.83 MeV level of 


2 ieee) Wy 
o as evidenced by the decided difference in character of the 
charged=particle spectra corresponding to the adjacent neutron 
time-of-flight channels (see Figure 3-17). The O14 group would 

be lost in the large peak at .75 MeV which is the result of random 
coincidences between neutrons and scattered deuterons from the 
beam. The hump in the spectrum at about 1.3 MeV could be a result 
of An and 45 summing in the detector. 

Curves a, b, c, and d in Figure 3-18 are unnormalized three- 
body energy distributions, one of which would be expected to 
obtain if the 11.83 MeV level disintegrates by a direct three- 
body process, Curve 'a' is the classical distribution for this 


case and curve 'b' is the classical distribution modified by 


penetrability factors for the three alpha particles. That is, 
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9 C2 a). 

f,(€1) = f (e))P, P P 
where the assumptions have been made that particle 1 has the angular 
momentum of the carbon state (J = 1), and that for an assigned energy 
Ey for particle 1, the other two particles share equally the remain- 
ing energy. Curve 'b' differs little in shape from 'a', both 
indicating that there would be an equal sharing of the disintegra- 
tion energy among the three alpha particles. Curves 'c' and 'd' 
are the Delves (loc, cit.) distributions with 2 


= 1, 2, = 0, and 


1 2 


dy = so) = l respectively, Explicitly, 


f(e) = eft - ¢yl/2 


and f , fe) = ell2y ~¢)°/? : 


Of these three-body distributions, only curve 'd' could be considered 
to account in part for the measured distribution, with the three 
peaks riding on it. However, if this were the case, the number of 
counts in the "random" peaks (caused by chance coincidences between 
neutrons and scattered deuterons from the beam) of the charged- 
particle spectra would not follow the neutron time-of-flight spectrum. 
Therefore, it is considered that the three-body distributions cannot 
account for the shape of the charged-particle spectrum, On this 
basis it is concluded that the 11.83 MeV state decays entirely 


through Be’, and further that, from a crude measurement of the areas 
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under the peaks at 2,75 (25) and 1.5 (A541 or 34) MeV, it is estimated 
that roughly ten percent of the transitions proceed via the oa 
ground state, 

The peaks in spectra no, 6 and 7 of Figure 3-17 at 2.25 MeV 
could only correspond to alpha particles from the 10,84 MeV state in 
cl? decaying to the ground state of Be?, This is in agreement with 


the only other observation of the decay of this state (see Chapter I). 


IV. Measurements Using Two Solid State Detectors 


A series of measurements employing two solid state detectors 
was carried out in an effort to detect two or all three of the alpha 
particles produced in the breakup of a given level of ae The 
counter geometry is shown in Figure 3-4b. The solid state detector 
(a) on the same side of the deuteron beam as the neutron spectrome- 
ter was used to provide the timing pulses, The pulse-height spectrum 
corresponding to the charged particles detected by the other solid 


state counter, a in coincidence with particles detected by a 


E’ T 
(with a resolving time of 35 nanoseconds), was analyzed with the 
neutron spectrum by the two-parameter analyzer set up as previously 
described, The complete electronic arrangement is shown in Figure 
2-5, Measurements were taken for several combinations of the 
detector positions, Ca.) ), at each of the neutron detector 
positions (0° and 45°), The running time for each geometry was the 


order of eight hours so that, in order to account for any drifts in 


the electronics during a run, the kicksorter was read out and cleared 


~ v2 - 


hejsemizas at sf .veM (re? xo ro”) Gof baw (ope) 8 I | | 
| Sa add siv bssc0rq emotyienst? sd3 to jnsorsiq: 


ye w ig 


Vem 25,8 os Vi-€ oxughd Yo T bee soe 


ft singe Vel BB101 ott most estotsing adate 0 
A9e dnoasdrgs nt st oiat ae thinset tnliid 
. (I yodqed9 see) otsta ality Yo raceb as Yo notsmrnande tase 


- a 


a 


etodo9i1sb sjsta bliiee ows oiiiaie ssnamrunai Yo 
eiigie stit Yo seria {Is 10 ows Joateb o3 d10%29 na. al a 


5) ofT, .~'9 lo [oval nevia « 0 quaesxd edo mb. 
wos2SI5b SBinse biloe siT .di=€ sruight ai svoda at ytemoy 
“smorjooqa mortusn si3 aa masd nors3ueb oi Yo abe sane 
switsaqe trgtet~safva aiT .eseluq gatats ort sbiverq o2 a 
| bhfoe tefjo. sd2 yd betost9b polotsuns, bagando ae 
| 7p? Xd bedseseb eslotitaq sat, sonabtsntos tage 
| ad) rittw besvians eaw “Cobnasseanea é | a 


ylavotvaerq 28 qv joe rssyigaa 3 FP es 


stug!t at nwore.et sroaegneny 2 
oe o- 


93 Yo arnotisntdmos £ 


tos2aI3eb nozjven edt te Bish 


at 


_ 


ott 2aw yriemeo3 fioss 103 oataian 
nl estizs yar 10% tmuonas oF 


} betse!o ban juo baat esw Fe 


=\re 


3S Ze 


TWO PARAMETER DATA FOR THE ALPHA 
PARTKLE BREAK-UP OF LEVELS OF 


1 
c* Feo By THE 8 '(4,n)c* REACTION 


@- 0 , Q@-90" @= 90° 
A ie M 


‘ | METER 


Figure 3-19, Isometric display of the two-parameter data obtained for 
@..2 09. So 8 1ometer Opie 900 Oe e000 Reale 


ee ts 


"OZ € e4nby 4 


}oOWosS | 


iA 


8 | 

) 

u 

U 

ul 
qHoit4 ~ 50 


me “oSt = U6 


4ed-om} ay} JO Ae, dsip = 


Me *sa}ow | 


“006 


Py ‘05g 2 406 
JO} paeulrezgo ejep sadqjowe 


“APW GS" I 


NEUTRON COUNTS (— 500) 
3 4 5 6 


———j———= v 
CHARGED PARTICLE 


- anil 


joes or fr 
e281 2 
Die eM2 8 


$.D 28 via 3° 
8 “Mvad Tome 1! 9 
u 


sD 
e-< 94335'"p wo @oezcen 0 


Gi-€ Oty 39S (Sp . 


@enu2aiy 


a) 
v 
m 
°o 
“4 
» 
= 
= 


WOLLVYRIZILN3O! ¥Y3d 


- 58 - 


after every 100 uc of charge collected at the target (i.e., about 
every twenty minutes). A check on the random counting rate was 


obtained by delaying the Ot pulses with respect to the a. pulses 


E 
by 70 nanoseconds, On the basis of 100 uc runs, the true to 
chance ratio was greater than 25:1 for every geometry. 

If the states of cl? decay sequentially through nei the 
kinematics are well defined, and peaks would be expected in the 
charged=-particle spectra. On the other hand, flat spectra would 
be expected if the levels disintegrate by a direct three-body 
process, The experimental spectra obtained for every geometry 
investigated exhibited considerable structure. Representative 


of the data are Figures 3-19 and 3-20. The former shows the 


data in isometric display for 6 = 0°, 6 = 90°, and 6 = 90°; 
n Oe Op 


and the latter for b = 45°, a = 85°, and as = 90°, The 
various peaks in these spectra can be attributed to alpha par- 
ticles originating from the 9,64, 10.84, 11.83, and 12.71 MeV 
states of of, and are labelled accordingly. The number 0 or 
2.9 in brackets following the cl? assignment indicates the state 
of He involved, These assignments were made on the basis of 
kinematic calculations performed by the IBM 1620 computer to 
determine the energies of the three alpha particles for each 
decay as functions of their laboratory angles of emission, 
assuming that the decay proceeds sequentially through Bar 


These calculations are fully described in Appendix I, Nuclear 


Reactions III. 
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V. Conclusions 


The alpha-particle energy spectra and angular distributions 
obtained in this work indicate that the 12.71 and 11.83 MeV levels 
of ake decay predominantly through ae The absence of transitions 
to the ited ground state supports the 1" spin-parity assignment 
for the 12,71 MeV state, About 90 per cent of the transitions from 
the 11.83 MeV state go through the 2.9 MeV state while the 
remaining 10 per cent feed the ground state. 

It is likely that the 10,84 MeV state of ct decays primarily 
to the ground state of Be’, This is in agreement with the only 
other observation of the decay of this state, Transitions to the 
first excited state of Be® are not ruled out by the present work, 
however, 

The angular distribution of the alpha particles from the 
9,64 MeV state strongly supports the 3 assignment for this level 
as opposed to the 1 assignment, 

The information obtained from this work on the alpha-particle 
decay of the levels between 9 and 13 MeV excitation, coupled 
together with the available information (see Chapter I) on the 
decay of the other levels in e supports the fe ag Two~-Body 
Cluster Model representation of the a nucleus, at least to 16 
MeV excitation, 

This work also eliminates the simultaneous three-body breakup 


of the carbon states between 9 and 13 MeV as being the source of 
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the neutron continuum observed in the neutron time-of-flight 
spectrum (see Figure 3-1), It is suggested then that this con- 
tinuum may be caused by the overlapping of the g levels including 
the very broad level at 10,1 MeV (Ni 62), Neutrons from the 
breakup of the 3.04 and the 2,43 MeV states of Be” formed by 


* 
the i tae ee” reaction will also contribute to this continuum, 
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APPENDIX I 
THE KINEMATICS OF NUCLEAR REACTIONS 


tT. Introduction 

Three programs have been written for the University of Alberta 
IBM 1620 computer to provide nonrelativistic solutions to the kinematics 
of nuclear reactions. The first and most generally useful program deals 


with reactions of the type 


M) + M, ny + my (1) 


or more briefly, Mo(My,M3)My. Here M, is the bombarding particle, My the 
target nucleus, M3 the light reaction product, and My the heavy (recoil) 

reaction product. The second program considers reactions of the type (1) 
where the residual nucleus, My, breaks up into two more particles, Me and 


M This program then is suitable only for cases when My, is formed in an 


6" 
unbound state. The third program goes one step further and determines the 
energies of the particles formed in the breakup of M¢ into M7 and Mg. This 
program was written specifically to provide the energies and angles of emis- 
sion of the three alpha particles produced in the breakup of os levels 


12* 
formed by the Ba sn)c reaction, assuming a two stage disintegration 


8 ’ 
process through states in Be , this reaction being the subject material of 


this thesis. 
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A. Nuclear Reaction I 

This program is concerned with reactions of the type (1). Ball (Ba62) 
has written a similar program for the 709/7090 computer at Oak Ridge, and 
some of the features of his treatment are incorporated here. Tables are 
also available (Ma59) for the computation of the energies of nuclear reac- 
tion products, and for the transformation of quantities between the labora- 
tory and the center of mass systems. However, tables are awkward and time 
consuming when complete kinematic solutions are required for a large number 
of reactions, as is the case when unravelling the complex spectra obtained 
in charged particle spectrometry. The program of Ball (Ba62) was unsatis- 
factory from the point of view of incompatibility with the 1620 computer. 

Nuclear Reactions I computes as a function of the laboratory angle of 
emission of the light reaction product, My» the corresponding center of 
mass angle, the laboratory energy of My the ratio of the cross-section in 
the center of mass to that in the laboratory for My» the energy and labora- 
tory angle of emission of the residual nucleus, Mis the momentum transfer 
wave number and the magnetic rigidity of My. An interaction radius parame- 


ter R/R, = mt/3 si wi/3 


l , ) is used to calculate the scattering argument. 


Input Data 


The input data arepunched on five IBM, 80 column, cards, The 
following table (Table 1) indicates what information is contained on each 
card, and at what position on the card the data should appear. The cards 
should be arranged in the order given below, Any number of sets of data 


may be stacked together. 
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Table 1 


Definition 


Title Card; the first 55 
columns, beginning with 
column 2, are available, 
Any notation to identify 
the reaction is permitted. 


Mass of the incident 
particle, My 


Mass of the target nucleus, 
Mo 

Mass of the light reaction 
product, M4 


Mass of the residual 
nucleus, M, 


The Q of the reaction to 
excited state S of the 


residual nucleus 


"Starting' incident energy 
of My 


Increment in incident energy 


Final incident energy of My 


Starting laboratory angle for 


particle My ( >0) 
Angular increment (>0) 


Maximum laboratory angle for 
M,(< 180°) 


Field 
begins at 
column 


amu 


amu 


amu 


amu 


Mev 


Mev 


Mev 


Mev 


radians 


radians 


radians 
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Output 


The listing of computed, angle dependent quantities is preceded 
by a heading which contains the input data (excluding the angle input 
data) and some computed, angle independent quantities. Table 2 indi- 


cates the signficance of the computed quantities in the heading. 


Table 2 
Symbol Definition 
The energy of the incident particle Mev 
(Mj) in the initial center of mass 
system 
EICOM The total kinetic energy in the initial Mev 


center of mass system 


EFCOM The total kinetic energy in the final Mev 
center of mass system 


E3COM The energy of M3 in the final center of Mev 
mass system 


E4COM The energy of M, in the final center of Mev 
mass system 


K1,K2 Wave numbers of My and M3 in the center | fermis 
of mass system 


R OVER RO The interaction radius divided by the 
nuclear radius parameter 


Table 3 shows the column headings for the listing of the remaining 


computed quantities. 
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Symbol 


THETALAB 
| THETACOM 

PHI 

E3 

E4 


XSEC-CMTOL 


KCOM 


KR OVER RO 


BRZ 


re. 
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Table 3 


Definition 


The laboratory angle theta of emission 
of M3 (input data) 


The center of mass angle corresponding 
to the input angle theta 


The laboratory angle of the recoil 
nucleus, My 


Laboratory energy of the light reac- 
tion product, M3 


Laboratory energy of the heavy reac~ 
tion product, My 


The ratio of the cross~section in the 
center of mass to that in the labora~ 
tory system for the light reaction 
product, M3 


The relative wave-vector of the par- 
ticles M; and M3 in the center of 
mass system 


The scattering argument divided by the 
nuclear radius parameter 


The magnetic rigidity of particle M3, 
the light reaction product 


Units 


degrees 


degrees 


degrees 


Mev 


Mev 


fermis 


inches 
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A typical Fortran II output page is shown in Figure 1 for the 


* 
reaction ee tae (12.71) at a deuteron bombarding energy, E,, of 


1.5 Mev. 


When the velocity of the center of mass in the initial system 
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exceeds the center of mass velocity of the light reaction product, Ms 
E3 becomes double valued, and the maximum allowable angle, theta, is 
restricted to values below 180 degrees. In such cases, the program will 
automatically replace THETF by the calculated maximum angle, provided, 
of course, that THETF exceeded this maximum angle in the first place. 
The two values of those items tabulated in Table 3 appear one after the 
other in the listing of the output, 

The computer time required for this program is three minutes per 


Q value per El value for theta varying from 0 to 180° in 5° intervals. 


Kinematics 
Figures 2 and 3 show the velocity diagrams of reactions of the 
type M, (M, »M,)M, « Figure 2(a) depicts the situation in the initial 


laboratory system with the bombarding particle M, approaching the sta~ 


1 


tionary target nucleus M, with a velocity v, and energy E Figure 


z 1 i: 


2(b) shows the reaction products in the final laboratory system in 
which the light product My, with velocity V3 and energy E, is emitted 
at an angle 6 with respect to the direction of the incident beam and in 
which the residual or recoil (heavy) nucleus, My? is emitted at an 
angle ¢ with velocity Y, and energy Eye 

Figure 3(a) shows the situation just prior to the reaction in the 
initial center of mass system. The projectile and the target nucleus 
are approaching one another with velocities v4 and Vo respectively, 


where 


Vv OD) 
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where ed is the laboratory velocity of the center of mass of the 


initial system, It is worth noting, for future reference, that 


i cha (4) 


E, (5) 


The total kinetic energy available in the initial center of mass 


system, E1COM, is given by 


2 grat 4 
EICOM = 1/2 MV, + 1/2 MV, = Mi, ri M, Ey (6) 
or more simply by 
M, 
EICOM = Ey se Be = M re M, Ey 


where Eis the energy of the center of mass of the initial system, 
cm 
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and is seen from equation (4) to be 


Es. «4 (7) 


The total disintegration energy in the final center of mass 
system, shown in Figure 3(b), is the Q of the reaction plus the 
‘inetic energy available in the initial center of mass system, That 


is, 


EFCOM = Q. + E1COM (8) 


With reference to Figure 3(b), and by applying the laws of con- 
servation of energy and momentum, the center of mass energies of the 


reaction products can be determined: 


2 2 

1/2 MV, + 1/2 M,V) = EFCOM (9) 

- 10 

and MV, = M,V)> (10) 
My 

from which E3COM = EFCOM Ci) 

M, + M, 

oe 

and E4COM = F—4 EFCOM , (12) 
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Convenient expressions for the laboratory energies of the reaction 
products, M, and Mis are obtained by applying the cosine law to the 


velocity triangles of Figure 3(b). The expression for E, is: 


z 2 2 
1/2 MV, 1/2 MV on + 1/2 MV, aa 2(M,/2)V__V,cos (0) (13) 
2 2 
or E, = A +B + 2A B cosO (14) 
3 fe) fe) 0 O 
M,M.E 
where a = pl 5 and B2 = E3COM, 
Cu, 7+") a 


In a similar manner, using the triangle of Figure 3(b) associated with 


residual nucleus, M)3 


2 2 
Ey, = Ay + By - 2A,B, cos (15) 
M.M,E 
where AS = Lo and BY = E4COM, 


Before equations (14) and (15) can be used, the input angle, 6, 
must be converted to the corresponding center of mass angle, 0, 
Applying the law of sines to the upper velocity triangle of Figure 


3(b), 


= = G 16) 
Sin(O - 6) = (V__/V,)Siné (16) 
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The symbol used in the program for the difference (0 -6) is D, and 
vo WN = Ao /B = y. Because the computer subroutine library does not 
contain the inverse sine functions, a transformation to a tangent 


function is necessary: 
Tan(D) = Sin(D)/{1 - sin?(p)}+/? (17) 


id, hence, © = ArcTan{Tan(D)} +6, (18) 


The angle of emission of the residual nucleus in the laboratory 
system can be obtained by applying the sine law to the lower velocity 


triangle of Figure 3(b). That is, 
Sing = (V/V oy) Sin (0 + o) 
or, by expanding the right-hand side and dividing by Cos, 


7 SinO 
A /By - Cos0 


(19) 


since Ven! V4 = A,/B,: 


Under certain conditions the energies of the reaction products 
will be double valued, and the laboratory angle, 6, will be restric~ 
ted to values less than 180°. These conditions are depicted in 
Figure 4, Figure 4(a) shows the usual case where the velocity of 
the center of mass, Mee is less than the velocity of Mas Va» in the 


final center of mass system, For this case, rE. is single valued, and 
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6 can take on any value up to and including 180°, When Vim tS equal 

to Vas as shown in Figure 4(b), the maximum value that 6 can take is 90° 
as this corresponds to 180° in the center of mass system, E, is still 
single valued, however, Figure 4(c) shows the situation where De is 
larger than V4! there is a second center of mass angle, On» correspon- 


ding to a given 8, and hence two values of E, for a given 6. The value 


3 
of 0, can be obtained from the geometry of Figure 4(c), and is given by 


e,=1+26-06,. (20) 
This new value for the center of mass angle is then used to recalculate 
all those quantities in the output listing which are angle dependent. 
The two values of each item follow one another consecutively in the 
listing. It can also be seen from Figure 4(c), or from equation (28), 
that a is restricted to a value less than 180°, and that it is 


given by 


Sine ox = V/V om = BA/AS = 1/y (21) 


from which 


1/2 


tS) = Rredddt tity - 1) 4 (22) 


max 


The computer examines the value of y to ascertain which of the above three 


conditions obtains: 
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j if y < 1, the computer proceeds with the range of 
as specified in the input data, 


ii, if y = 1, THETF is replaced by 1/2 provided that 
THETF exceeded 1/2 in the first place, 


» Ss Se © ee ee Base is computed from equation (21), 
and replaces THETF provided that THETF exceeded 


s) in the first place, 
max 


The Center of Mass Correction for My 


Before measured cross sections or angular distributions can 
be compared with theory, it is necessary to convert them to the center 
of mass system of coordinates. The total intensity of reaction 


product, tas is the same in both systems, so that 


01 ap 60) Sin6de = 6,,,(0)SinOdd (23) 


CM 


where 01 op 69) is the differential cross section as measured in the 
laboratory system, and Foy 69) is the differential cross section as 
measured in the center of mass system. The azimuthal angle is the 


same in both systems and therefore has not been included in equation 


(23). The correction factor required is 


Soy) Sing do 5 
= es 2 
O1 ap 69) SinO 
With reference to Figure 3(b) 
Sint = Vi/v (25) 
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and Sin(O - 6) = (V om/ V3) Sine = ySiné (26) 
from which Se a ~ {sin™ (ySiné)} 


2 2 e172 
so that dé (1 aie a! ) 
~ n 6 (27) 


ag yCos6 + (1 - y"sin-g) Lf? 


However, by applying the cosine law to the upper velocity triangle 
of Figure 3(b) it can be seen that the denominator of equation (27) is 
simply V3/V 3. That is, 


/2 


v, = V,[yCosé +(1 - +fin ey" 1 « (28) 


Therefore, Sadie? a ae y*sin29)1/? 


97 ap 69) v2 /Vv2 


= (B°/E,)(1 - y*sin’a)}/? (29) 
Equation (29) is in terms of previously computed quantities, and is 
the one which is used to calculate the center of mass correction 
factor, 
Other quantities which are useful in angular distribution analysis 
are the wave numbers of the incoming projectile and the light reaction 


product, M,, and the relative wave number of the transferred momentum. 


3 


These can be obtained directly from quantities which have already been 


computed, The general expression for wave number of a particle of 
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mass M and with energy E is 


1/2 
K = ss faynt (30) 


where fi is Planck's constant. Using this relationship, the equation 


for the wave number of My in the initial center of mass system is 


2 Z 


K} = M,(E1COM)/20.908 fermi ~ . (31) 
And for M, in the final center of mass system, 
2 ~-2 
Ky = M, (E3COM) /20 .908 fermi 0 (32) 


The relative wave number, K, of the transferred momentum is then 


given as shown in Figure 5 by 


2 2 Z 
K = K + Ky = 2K K,Cos0 (33) 


so that KCOM = K feet. 


Figure 5: Wave vector diagram, 
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The final parameter listed in the output is the magnetic rigidity 


of the light reaction product, My. This is given by 


BRZ = 56.67(MgEq) /7 kilogauss-inches (34) 


where B is the magnetic field in kilogauss, R the radius of curvature 


in inches, and Z the electronic charge of M It was felt that this 


3° 
parameter would be very useful in connection with the broad range 
spectrograph and spectrometer facility being set up at the University 
of Alberta, 

The computer program for the foregoing was written in FORGO 


and then compiled in FORTRAN II, The FORGO listing is presented on 


pages /6 and 77, 
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ee IDENT 29031659KINETICS OF NUCLEAR REACTIONS »WeCeOLSEN 


C 


100 
1 


102 


103 
104 


105 
106 
107 


row £ 


We CeOLSENs903165sKINETICS OF NUCLEAR REACTIONS I 
FORMAT (55H 
REAN 100 

PUNCH 100 

READ 1019EM1sEM2 9EM3sEM4 

FORMAT(1X4F12.6) 

READ 101%:QS 

PUNCH 102%EM1sEM29EM39EM4 

FORMAT (4HOM1=F10.6 3X9 3HM2=F 10069 3X9 3HM3=F 10069 3X9 3HM4=F 1006) 
PUNCH 1109QS 

FORMAT (4HOQS=F10.6) 

READ 101sE1S9DE1 sE1F 

READ 101sTHETSsDTHET»sTHETF 

El=E£1S 

FI1COM=E1*(EM2/(EM1+EM2) )**2 
EICOM=E1*(EM2/(EM1+EM2) ) 

EFCOM=QS+EICOM 

AO2=EM1*EM3#E1/( EM1+EM2) **2 

BO2=(EM4/ (EM3+EM4) )*EFCOM 

AO=SQRTF(AO2) 

BO=SQRTF(BO2) 

AL2=EM1L*¥EM4*E1/(EM1+EM2) **2 
B12=(EM3/(EM3+EM4) ) *EFCOM 

Al=SQRTF(A12) 

Bl=SQRTF(B12) 

ROVRO=(EM1¥* 9 333334+EM2** « 33333) 
CAY12=EM1*E1COM/20.908 

CAY1=SQRTF(CAY12) 

CAY322=EM3*B02/20.908 

CAY3=SQRTFI(CAY32) 

PUNCH 103sE1sE1COMsEICOMsEFCOM 

PUNCH 1049B8B029B12 

FORMAT (GHOE1L=F9 059 3X96HELCOM=F 9059 3X 9 GHE I COM=F9 659 3X9 6HEFCOM=F 965 ) 
FORMAT ( THOE3COM=F965 9 3X9 6HE4SCOM=F 965 ) 
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THETFH=ATANF(C/SQRTF(1e¢-C¥*¥2)) 

GO TO 5 
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THETA=THETS 
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SIND=GAM*SINF (THETA) 

TAND=SIND/SORTF(1e-SIND¥*2) 
COMRD=ATANF(TAND)+THETA 

COMTH=COMRD*¥5722957786 

F3=A024+R02+26¢*AO*BO*COSF (COMRD ) 
F4=A12+B12—-2+*A1*B1 *COSF (COMRD ) 

TANPH=SINF (COMRD)/(A1/B1-COSF (COMRD) } 
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PHIRD=ATANF ( TANPH) 

IF (TANPH) 83899 

PHIRD=3e1415927+PHIRD 
PHI=PHIRD¥5722957786 

GAM2=GAM**2 

SNTH2=SINF (THETA) *¥¥*2 
SIGR=BO02*SQRTF(1e¢-GAM2*SNTH2)/E3 
CAY2=CAY12+CAY32-2 -*CAY1¥CAY 3*¥COSF (COMRD ) 
CAY#SGRIFICAY2) 

CAYR=CAY*ROVRO 

BRZ=56e67*SQRITF (EM3*E3) 

PUNCH 108s THETL»COMTHsPHI sE39E4sSIGRsCAY sCAYRs BRZ 
FORMAT(FIO0e292Flle292Flle3sFl2e4/1HI2Z4XFBe49F1l2049F 1404) 
TF (C#1.)109813913 

COMRD=3-1415927+2.* THETA-COMRD 
COMTH=COMRD*5722957786 
E3=A02+B02+2¢*AO*BO*¥COSF (COMRD ) 
E4=A12+B12—-2¢%*A1*B1*COSF (COMRD ) 
TANPH=SINF(COMRD)/(A1/B1-COSF (COMRD) ) 
PHIRD=ATANF(TANPH) 

IF (TANPH) 11911912 

PHIRD=321415927+PHIRD 
PHI=PHIRD*5722957786 

SIGR=BO02*SQRTF (1+¢-GAM2*SNTH2)/E3 
CAY2=CAY12+CAY 32-2 e ¥CAY1L*¥CAY 3*COSF (COMRD) 
CAY=SORTFI(CAY2) 

CAYR=CAY*ROVRO 

BRZ=56-267*SQRTF(EM3*E3) 

PUNCH 108s THETL»COMTHsPHI 9E39E49SIGR sCAY sCAYRsBRZ 
THE TA=THETA+DTHET 

IF (THETA-THETF)696914 

E1=E1+DEl1 

FETEL=E1F 315915916 

PUNCH 109 

FORMAT(1H2) 

GO TO 2 

PUNCH 109 

GO TO l 

END 
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B, Nuclear Reactions II 


This program is an extension of Nuclear Reactions I in the 
sense that it considers those reactions in which the residual nucleus. 
Mis is left in an unbound state which decays, in the reaction plane 


defined by M, and Mie into two more particles and/or nuclei, Me and 


Me? The coplanar velocity diagram of the final system is shown in 
Figure 6, The program instructs the computer to determine the labo- 


ratory energies and angles of emission, (E, yy) and (E,.B), of M. and 


Me as well as Ey, 


Reactions II also determines, for Mey the angles of emission, yp, in 


Oo, Ey, and ¢ defined in Nuclear Reactions I. Nuclear 


the center of mass system of M,, corresponding to the angles y, and 
the ratio of the cross section in the center of mass of My to that 


in the laboratory system, 


Input Data 


Seven standard, 80 column IBM cards are required to accommodate 
one set of input data for this program, Table 4 indicates the 
contents of each card. Any number of sets of data may be stacked 


together, each separated by cards 8 and 9 of Table 4, 


Output 


A typical output page is shown in Figure 7a for the reaction 
(A nC (9464) >art Be°(0). The listing is self-explanatory 


except perhaps for the following symbols: 
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B11(D»N)C12(9.64) TO ALPHA PLUS BE8(0) 
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Table 4 


_oymbol . Definition 


This is the Title Card, 
The first 55 columns, 
beginning at column 2, are 
available, Any notation 
to identify the reaction 
is permitted, 


The mass of the bombarding 
particle, My 


The mass of the target 
nucleus, M, 


The mass of the primary, 
light reaction product, My 
The mass of the primary 


residual nucleus, My, 


The mass of the secondary, 
light reaction product, M. 


The mass of the secondary, 
heavy reaction product, Me 


The Q of the primary reac~ 
tion to the excited state, 


So, Or fhe My nucleus 


The disintegration energy 
of the secondary reaction 
in which M, breaks up into 
Me and Me 

The starting, or initial 
laboratory angle of 
emission, 6, of Me 


THETS 


DTHET The increment in 6 


THETF 


The final or maximum angle, 
6, of interest 


At 


Begins — 


Column 


14 


26 


ho 


14 
26 


Units 
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amu 


amu 
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Table 4 (cont'd) 


Begins 
At 

Column Units 
The starting angle of emis- rad, 
sion, p, of M, measured 
with respect fo the labora- 
tory direction of the 
residual nucleus, M,, as 


shown in Figure 5, This 
is usually zero, 


The increment in the angle rad, 
D » 


The final or largest value rad, 
of p of interest. This is 
usually 27, 


The initial value of the Mev 
bombarding energy of M) 

The increment in El Mev 
The final value of the bom~ Mev 
barding energy of My 


eol 


. sBRANCH O040R 


i, THETA-COM,. This is the center of mass angle, 0, 


corresponding to 6, 


ii, XSEC-CMTOL. This is the notation used for the ratio 
of the cross section for Me in the center 
of mass system of M), to that in the 


laboratory system, 


All mass values are in amu, all angles in degrees, and all energies 


in 
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Kinematics 


The kinetic energies of the products of the secondary reaction 


are, in the center of mass system of Mie 


6 
5 = Oi +t) Qn , for Me 5 (34) 


M 


where Q, is the disintegration energy of the secondary process, 


Expressions for the laboratory energies of Me and Me can be 
derived with the help of Figure 6, Applying the law of cosines to 
the velocity triangle bounded by Vy» Vos and Ve one gets 


2 2 Z 
V5 a ey + Ve > 2v,V.Cosp ‘ (36) 


and upon multiplying through by M,/2 


/2 


1 
= : = Ee FeE 37 
Ee ESF + E.P 2(E.F E.P) Cosp (37) 
where Vv, = the laboratory velocity of Me 
Ve = the velocity of Me in the center of mass 
system of M,, 
iy = } 5 . 
ESF (Mo /M) JE, 


E, is given by equation (14), 


Considering the triangle bounded by the velocities Vi» Veo and 
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Ve of Figure 6 and again applying the cosine law, one gets 
2 2 


2 
Ve = Vv, + Ve + 2v ,V.Cosp (38) 


and by multiplying through by Me /2 


= Lé2 
E EGF oe EP + 2(E (FE (P) Coso (39) 
where Ye = the laboratory velocity of Me 
= the velocity of Me in the center of mass 


system of M 


EGF = (Me /M) DE, 


A 


Equations (37) and (39) are used to calculate Es and Ey as 
functions of the angle, 0, between Vy, and Ve» 
It is worth noting here that op is an arbitrary angle which forms 


part of the input data, and which varies from zero to 360° in small 


increments of 10 or 15°, say. 
The angles of emission of Me and Me can be determined by first 
calculating the values of 6 and ¢« for each value of p. A convenient 


expression for 6 in terms of previously determined quantities is 


Gong iw eee TES <P (40) 


For computer convenience 6 is denoted by DL in the program, 


For p < m, the laboratory emission angle of Me, y, is given by 


Y= ¢- 6 (41) 
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¥ * Ge Oe (42) 
Epsilon, the angle between Vy, and Veo is given by 


(E, + E,F - E,P) 


Cose «= , (43) 
2(E, EF)! 
6 6 
For p < 7, the laboratory emission angle of Mes 8B, is given by 
B= o+eE (44) 
and for p > 7, B= ¢-e, (45) 


The Center of Mass Correction for Me 


Consider the laboratory system of coordinates (x,y,z) of 
velocity space as having its z-axis in the laboratory direction of 
residual nucleus, Mae of the primary reaction. That is, the z-axis 
makes an angle @ with the direction of the incident beam of bombar- 
ding particles. And let the coordinate system (x', y'’, z') in 


which M,, is at rest have its origin at the end of the velocity 


vector ae and its z'=axis in the direction of the velocity of 
M) in the center of mass system of the primary reaction (i.e., V,)« 


These two coordinate systems are shown in Figure 6, This choice of 


> 


the z'=-axis is convenient when considering the angular distribution 


of particles from the breakup of levels formed in a stripping 
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reaction, for example, the oe we reaction. If pure Butler 
theory applies, the angular distribution of the a-particles from the 
decay of ee leading to states in Be® will be symmetric about the 
direction of the captured proton; that is, about the z'-axis, 

The primed system is related to the unprimed system by a 
translation, Ugo and a rotation of a = (nm - 0 -¢ ), both in the 
(y,z) plane, the plane of the primary reaction, If (6,n) and (y,v) 
are the polar and azimuthal angles in the unprimed and primed 


systems respectively, the transformation equations are 


V,SinyCosv = v,SinéCosn (46a) 
V.SinySiny = v,SinéSinnCosa - (v.Cos6 - v,,) Sina (46b) 
V.Cos = v,SinéSinnSina + (v,Coséd - Vv) Cosa. (46c) 


The total number of particles, Mes detected is the same in both 


systems. Therefore, 
0(6)Sinédédn = o(p)Sinypdydy (47) 


where o(6) and o(w) are the cross sections for the production of Me 
in the laboratory (unprimed) system and the center of mass system 


(primed) of M, respectively. The required correction factor is given 
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by 
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Now, dé = BV 505 + aye! + dv (49) 
_ 23 3 3 
and dn = ave %5 rs Sydv of Slay ; (50) 


Ve is a constant quantity and therefore dV. = 0, The partial deriva- 


5 
tives 36/dv and dn/dv can be obtained by differentiating equations 
(46) partially with respect to v and solving the resulting system of 
three equations for the derivatives, In a similar manner, by 
differentiating with respect to wy, 96/aW and dn/dyp may be determined, 
The results of these calculations for Me confined to the primary 


reaction plane where n = v = #tn/2 are: 


06 on 
come = nce = l 
v ow Ug (51) 
V 
23s = cos(y - 6 #0) , (52) 
vy) Vv 
5 
V5 siny 
and On Lf md on ; (53) 
dv Sind 
5 
In equation (52), the plus sign corresponds ton = v = + n/2 
and the minus sign ton = v = -m1/2. 


The cross section ratio, equation (48), can now be rewritten as 
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From Figure 6 it can be seen that in terms of the arbitrary 


angle p, and for n = v # m/2, is given by 
9) a i eae (55) 


Forn = v ®& e=n/2, Vinj2 = ~Yqsg* In the computer program, v 
is computed from equation (55) only. Hence positive and negative 
values of p» correspond to cases when n and vw are +n/2 and -1/2 
respectively. Absolute value signs around \ in equation (54b) are 
required to eliminate this distinction. 

From the point of view of programming convenience, it is 
desirable to have an expression for the cross section ratio in 
terms of the arbitrary angle p. Such an expression can easily be 


obtained by noting that 


26 36 ap a6 
ay 0 BW ; see) 
since 39/ay = 1 from equation (55), and that, from Figure 6, 
ty 
Sind = -|=—=|Sino . (57) 
Vv 
5 
are Fe 
Hence 6 = ArcSin{-(E,P/E,) Sino} (58) 


Differentiating this expression with respect to p, one obtains 
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Using equation (37) for E 


59 
avE. VE PEF 
77 ee - Sinp (60) 
5 
E.P 1/2 SEF Sint = Cosp 
Hence» 28. = 2s a > 2 5 5 
dp p E,P > (E 3/2 
l- Eoin 0 5 
5 


(61) 


Substitution of equation (61) for 36/39 and equation (53) for 


dn/dv at n = v #& tm/2 into equation (54a) gives 


2 
oy) E.P | YE,P+E,F Sin’p - E,Cosp| is 
0 (6) 2 re FP: 
(E,) [1 - (E,P/E,)Sin’o] 


Equation (62) is used to calculate the cross section ratio for 
each value of the angle op, The results are tabulated under the 
column heading XSEC-CMTOL, The absolute value of the binomial term 
is taken since negative values of this ratio are unphysical, 

The program was written in FORGO and compiled in FORTRAN II, 
Five minutes of computer time are required to produce the output for 
one bombarding energy value (E,) and one angle 6, for o varying from 
0 to 2n in 10° intervals. The listing of the program statements is 


given on pages 88 and 89. 
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ee IDENT9903157sKINETICS OF NUCLEAR REACTIONSsTWOs WeOLSEN 
eeLOAD FORTRAN DUMP EXECUTE 


C 


100 


101 


107 
108 


WeOLSEN99031579KINETICS OF NUCLEAR RFACTIONSs ITI 
FORMAT(55H 

READ 100 

PUNCH 100 

READ 19019EM12EM2 sEM359EM49EM5 EMG 
FORMAT(1X3F12.6) 

READ 10159QS3QB 

READ 101lsTHETSsDTHETsTHETF 

READ 101sROST  sDROsROF 

READ 1O1sE1SsDE1+2E1F 

PUNCH 1029EM1sEM29EM3 

PUNCH 1039EM49EM5 5EM6 

PUNCH 1043QS5:98 

FORMAT (4HOM1=F10.693X9 3HM2=F10 069 3Xs 3HM3=F1006) 
FORMAT (GHOM4=F10.693Xs 3HMS5=F10669 3X9 3HM6=F1006) 
FORMAT (4HOQS=F10e623X»s 3HQB=F10-6) 

El=E1S 

PUNCH 105%9E1 

FORMAT (4HOE1=F10-¢6) 

AO2=EM1 *EM3*E1/ ( EM1L+EM2) ¥¥*2 

BO2=(EM4/ (EM3+EM4) )* (QS+EM2¥*E1/(EM1+EM2) ) 
AO=SQRTF(AO2) 

BO=SQRTF(BO2) 

THETA=THETS 

SIND=AO*SINF (THETA) /BO 
TAND=SIND/SORTF(1e~-SIND*¥*2) 
CAPTH=ATANF(CTAND)+THETA 
F3=A024+BO02+2.6*AO*BO*COSF (CAPTH) 
A12=EM1*EM4*E1/(EM1L+EM2 ) *¥2 

B12=(EM3/ (EM3+EM4) ) *(QS+EM2*E1/ (EM1+EM2) ) 
Al=SQRTF(A12) 

B1=SQRTF(B12) 

F4=A12+B12-2¢*A1*B1*COSF (CAPTH) 
TANPH=SINF(CAPTH)/(A1/B1-COSF (CAPTH) ) 
PHIRD=ATANF ( TANPH) 

IF (TANPH) 29293 

PHIRD=PHIRD+3.21415927 

PHI =PHIRD¥*5722957786 

THE TD=THETA*5722957786 

COMTH=CAPTH¥5722957786 

PUNCH 106sTHETD»COMTHsE3 

PUNCH 107sPHI5E4 

FORMAT ( 7HOTHETA=F 7029 3X9 LOHTHETA~COM=F 7029 3X9 3HE3=F 70 3) 
FORMAT (SHOPHI=F T7029 3X9 3HE4S=F 70 3//) 

FORMAT (4X 5HGAMMA6X3HPS I 8X2HE54X1OHXSEC-CMTOL4X4HBETATX2HES6 }) 


PUNCH 108 

ESP=EM6*QB/ (EM5+EM6 ) 
E6P=EM5*QB/(EM5+EM6) 

RO=ROST 

E5F=EM5*E4/EM4 
E5=F5F+ES5P—-2e*SQRTF(ESF*E5P) *COSF (RO) 
COSDL=(E5+E5F-E5P)/(26*SQRTFCESFH#ES) ) 
SNDOL2=(12¢-COSDL**2) 

IF (SNDL2)89999 
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SNDL2=0200000 
TANDL=SQRTF(SNDL2)/COSDL 

DEL TA=ATANF(TANDL) 

IF (TANDL) 10910911 

DEL TA=DELTA4+3.21415927 

IFARO- 261A 1 S927) 12812913 
GAMRD=PHIRD~DELTA 

GO TO 14 

GAMRD=PHIRD+DELTA 
GAMMA=GAMRD¥*5722957786 
E6F=EM6*E4/EM4 
E6=E6F+E6P+2e¢*SQRTIF(E6F*XE6BP) *COSF (RO) 
COSEP=(E6+E6F~E6P)/(2-*SQRTF(E6#E6F) ) 
SNEP2=(1¢-COSEP*%*2 ) 

IF (SNEP2)15916916 

SNEP2=0.00000 
TANEP=SOQRTFCSNEP2) /COSEP 
EP=ATANF(TANEP) 

IF (TANEP) Ts17918 

EP=EP+3.21415927 

IF (RO—3.61415927)19919520 
BETRD=PHIRD+EP 

GO. 10 21 

BETRD=PHIRD-EP 
BETA=BETRD¥57.2957786 
PSIRD=RO+CAPTH+PHIRD-622831854 
PSI=PSIRD*57.295/7786 

SINRO=SINF (RO) 

COSRO=COSF (RO) 

ES5R=ESP/E5 

P=E5R/E5 
Q@=}]./(SQRTF(1.-E5R*SINRO*¥*2 ) } 
R=ABSF(SQRTF(ES5SP#E5F )*SINRO#*2~E5*COSRO) 
SIGR=P#Q*R 

PUNCH 109sGAMMAsPSIsE5sSIGRsBETAsE6 
FORMAT(2F10.29F1l0e39Flle32F10e29F 1063) 
RO=RO+DRO 

IF (RO—ROF 149496 

PUNCH 110 

FORMAT (1H2) 

THETA=THETA+DTHET 

IF (THETA-THETF) 121922 

E1=E1+0E1 

IF(E1-E1F) 797923 

STOP 

END 
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Velocity diagram for Nuclear Reactions I11. 
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C. Nuclear Reactions III 


Nuclear Reactions III is again a straightforward extension of 
Nuclear Reactions II, In this case, Me is a nucleus formed in an 
unbound state by the first two stages of the reaction. It subse- 


quently breaks up into the two particles M_ and M,. It is the 


7 8 

purpose of this program to instruct the computer to determine the 
following quantities for a given value of E, and 0: (E,,0), (E, 29) 
(EosY)s Eos (E,,8,) and (E,58,)+ (E,,8,) and (Eos8,) are the la- 
boratory energies and angles of emission of the particles Mo and 
Mg respectively, The other variables have been defined in the 


previous programs, and all are shown in the velocity diagram of 


Figure 8, 


Input Data 


Six IBM cards comprise a set of input data for this program. 
Table 5 explains what information each card should contain, The 
ordering of the cards follows that indicated in the Table, The 


program must be read in for each set of data, 


Output 
A typical FORTRAN II output page is shown in Figure 9, The 
reaction in this example is: 
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Figure 9. Sample of Kinematics for the reaction: 
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Table 5 


Definition 


EM1 
EM2 
EM3 
EM4 


EM5 
EM6 
EM7 
EM8 


ROST 


DRO 
ROF 


ZETST 


THETA 


Mass of 


M 
Mass of M 
Mass of M 

M 


Mass of 


Mass of M 
Mass of 
Mass of 


Mass of M 


The initial angle of emission 
po of M. with respect to the 
laboratory direction of M,. 
See Figure 8, This is usuall 
zero, 


The increment in 0 


The final angle, op; usually 
27, 


The initial angle of emission, 
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for 6 = 45° and Ey = 1,5 Mev, 


The first two lines consist of the eight mass values, four per 
line, and the third line consists of the input energies Ey Qo» Qn» 
and Qn These quantities are in fixed point format. The next line 


records, from left to right; 6, 0, Eas ¢, and E The listing which 


4° 
follows consists of two parts: a line containing the values of 


vig Ee» and Ee for each value of 0, followed by four columns of 


single spaced numbers which correspond to, from left to right; Bie 


Es Bos and Eas These last four quantities correspond to ¢ as it is 


varied from its minimum to maximum values (0 t» 7 radians in 15° 
intervals in the above example). This format for the listing is 
repeated for each value of p as it is varied from its minimum to 
maximum values (6 to 27 in 15° intervals in the above example). 
Because of the periodicity of the trigonometric functions, By and 


Bo are correct to within +t nn, where n is an integer, Specifically, 


in the above example, 
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neutron 
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and considering the first row only in the listing for Mo and Mes 


ail S972" fern? =: 7120, 28% 9 ok = 0.139 Mev 
a2 


and B, = 419,779 = 27 = 59,77°, a 4,297 Mev. 
3 


Kinematics 


Only the kinematics for the breakup of Me need be considered 
here as the processes leading to the existence of Me have been 
treated previously, 


The kinetic energies of M_ and M, in the center of mass of 


v/ 8 
Me are 
Mp 
Ef 3 eee (64) 
Fi (Ma + Mg) r 
for Mo, and 
Mo 
fr UD ey ie 


for Mos where Qn is the disintegration energy of Mes 
An expression for E, in terms of the arbitrary (input) angle 
cS can be obtained by considering the velocity triangle bounded by 


Ve» Vos and V5 (Figure 8). Thus, 


2 2 2% = 66 
v7 = % + Vo 2v,V Cost ; (66) 


Multiplying through by (M,/2) 


1/2 
E, = EF + E.P - 2(EjF+E,P) lecosz ; (67) 
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For programming convenience 
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by velocity vectors Ve» Ves and Vo: Thus, 
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APPENDIX II 


A TRANSISTORIZED, FREQUENCY MODULATED PROTON MAGNETIC RESONANCE 
SYSTEM FOR THE MEASUREMENT AND CONTROL OF MAGNETIC FIELDS 


I, Introduction 


The original research project of the author was the instrumentation 
and utilization of a Spectromagnetic Industries 50 centimeter combined 
broad range magnetic spectrograph and spectrometer. At that time 
it was felt that the precise magnetic field regulation which can be 
obtained using nuclear magnetic resonance techniques would be 
necessary for the successful operation of this instrument. However, 
although there were several systems using this phenomenon for 
field regulation commercially available, not one had probes suitable 
for use with our magnet, and the purchase of custom, or specially 
designed probes, would have been unjustifiably costly. It was 
therefore decided to build our own regulating system, tailoring 
it to satisfy our particular needs, 

This NMR (nuclear magnetic resonance) system was designed to 
measure and control magnetic fields from one to thirteen kilogauss 
to a very high degree of accuracy. It is fully transistorized and 
employs frequency modulation of the radio frequency oscillator 


rather than the usual amplitude modulation of the magnetic field. 
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The circuits are of modular construction to facilitate trouble shooting 
and alignment procedures. The use of frequency modulation permits 
considerable reduction in probe size over that of conventional systems 
because the field modulating Helmholz coils are not required. This, 
in fact, was the main reason why frequency modulation was employed, 
It also gives the probe cylindrical symmetry so that special orienta- 
tion of it in the magnetic field to achieve the optimum signal to 
noise figure is not necessary. This feature has proved to be a great 
convenience, 

Four probes are required to cover the frequency range from 
4.8 to 55 mc/s ( 4 1 to 13 kilogauss), and they can be quickly and 
conveniently changed without breaking the vacuum of the spectrometer, 
The magnetic field is controlled by an error signal, derived from 
the field itself, which is used for closed loop regulation of the 
Magnet current power supply. It provides an error signal of +5 
volts for a 0.05 gauss change in magnetic field. } 

It should be mentioned that the homogeniety of the field 
produced by the spectromagnet was found, using the above system as 
a fluxmeter, to be far below contract specifications. The magnet 
was therefore returned to the manufacturer for rebuilding. It was 


for this reason that the author discontinued work on this project. 


Lbe Background and Design Considerations 


The operation of this instrument is based on the phenomenon of 
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nuclear magnetic resonance absorption (B148) in which radio frequency 
energy is transferred from an electrical circuit to a nuclear spin 
system located in a magnetic field. 

Consider an assembly of non-interacting nuclei in the presence 
of a magnetic field, Hoe For each nucleus there will be (2I + 1) 
equally spaced energy levels, where Ifi is the maximum z component 


of angular momentum, The energy of the ‘lis level is given by 


E,=E -(@.ii) 1 

ic Sey (1) 
> > > 

where yp is the magnetic moment of the nucleus, -yp . Hy the magnetic 

dipole interaction energy, and EG the energy of the unperturbed 

nucleus, Taking the z axis in the direction of the magnetic field, 


Hes E, may be rewritten as 


E,=E -u.Hk=E - tim (2) 


i fe) fe) re) 


since vi = yh a where y is the nuclear gyromagnetic ratio, The 


> > 
scaler product I .k is simply m the magnetic quantum number 


{? 
whose value runs from -I to +I in integral steps. Adjacent levels 


are thus separated in energy AE where 
AE = yhH . (>) 


This is the amount of energy that is required to induce transitions 
from a lower magnetic substate to an upper magnetic substate. It can 


be supplied to the spin system by placing the assembly of nuclei inside 
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a radio frequency coil which is driven at a frequency given by 


hv = 
v Y tH 
or 2mv = yH 
v= yH, (4) 
or = 
W YH, 


The nuclei used in the instrument to be described are hydrogen, or 


protons, for which 


Th = BS ~ 4.95784 mc/s/kilogauss (5) 
2m 4 mc 
where g = 5.58550 nuclear magnetons for protons 


e = electronic charge 
m_ = proton mass 


velocity of light 


QD 
" 


Therefore, equation (4) can be written, for protons, 
v = 4.25784 Hy Mc/s (6) 
with A, measure in kilogauss,. 


The accuracy to which the field can be measured, then, is limited only 
by the accuracy to which the absorption or resonance frequency can be 
measured, The frequency range of the rf oscillator is specified by 
this equation (6). Therefore, in order to measure fields in the one 
to thirteen kilogauss range, an oscillator with a frequency range from 


4,2 to approximately 55 Mc/s is required, 
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The optimum rf magnetic field that the oscillator should generate 
inside the coil can be estimated from the "saturation condition" 
(B148). That is 


y ae 4 
Y HjT,1,* = 1 (7) 


where Hy is one half the amplitude of the rf field, which may be 


written 


-2nivt 
H, = 2H, e » HL = 0 (8) 


T, is called the spin-lattice relaxation time, It is the time 
required for the spin temperature of the nuclei to return to the 
lattice temperature after the sudden removal of the rf field or, in 
other words, the time required for the state populations to return 


to their thermal equilibrium values by the interaction of the 


range from 


2 


spins with the lattice. The experimental values of T 


a to 104 seconds, being largest for solids, and as low as 10 


i 


seconds for pure liquids. The lowest values can be reached by the 
addition of paramagnetic ions, as will be mentioned later. Two 
very convenient proton sources, water and glycerin, have spin- 
lattice relaxation times of 2.3 and 0.007 seconds respectively. 
T,* is the characteristic time which describes the spin-spin 


interaction between neighboring dipoles, and is a direct measure of 


the inverse line width. It is defined by the equation 


T* = 1/2 8(Y) nag, (9) 
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where 1) Bie is the peak value of the observed shape, g(v), of the 
absorption line, Experimentally, g(v) is intermediate between a 


damped oscillator and a gaussian for aqueous ferric nitrate solutions. 


For a gaussian curve 


iy2 


1/T,* = (2m) yAH, (10) 
and for a damped oscillator 
1/T,* = (73/2) yAH (11) 


where AH is the width at half maximum of G(v), in gauss, 
An order of magnitude estimate of the natural line width may be 
obtained by considering that in addition to the main magnetic field, 


Ho each nucleus is subjected to 'local' field, H » set up by neigh- 


loc 


boring magnetic dipoles, In order of magnitude 
: = 1k 
H = u/r” = a few gauss for r = 1 (12) 


so that the resonance condition now becomes 


v= y/2n(H | a: ee (13) 


Loe 


The width of the absorption line in terms of frequency spread would 
2 a 7 g 
then be of the order of u /ir. or 10 kc/s for r= 1A, 
It should be noted that since T,* is an experimentally determined 
quantity, it necessarily contains all sources of line broadening, 


including inhomogenieties in the main magnetic field, Hoe 
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Bloembergen, Purcell and Pound (B148) give an expression, based on 
equation (7), for the optimum signal to noise ratio and it should be 
examined before using equation (7) to determine the optimum value of 


the rf field. This expression is 


2f3 222 
‘ V tYh Noi + 1) 


AJA eo 1462010 ) eee aes : 


14 1/2 
Y (TMT. (14) 
kT (kTBF) //2 ° 201 


where LP = the volume inside the rf coil, 


¢ = the fraction of that is occupied by the sample, 

Y = a constant, the order of unity, and depends on the exact 
line shape and the amplitude and shape of the modulating 
wave, 

No = the number of sample nuclei/cc present, 


= Boltzmann's constant, 


and T = the absolute temperature of the sample, 


B and F are parameters characteristic of the detecting apparatus, B being 
the band width of the amplifier system used, and F being an over-all 
noise figure. 

With equation (14) in mind, the proton sample chosen for the 
instrument to be described was water with ferric nitrate (Fe (NO) ,*9H,0) 
added in the proportion of one part to 100 parts of water. The 
addition of this paramagnetic salt reduced the value of T, from 2,3 
seconds to approximately 2 x 10° seconds, as determined from an 
experimental curve given in (B148). Extrapolation of figure 16 in (B148) 
gives the same value for T,* > Thus, the signal to noise ratio was 


increased by a factor of approximately 400 as compared to that for pure 
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water. Also, Ty) now forms an important contribution to the line width. 
Using equation (7) and T, = “= a2 & 10°? seconds as determined 


above, the optimum rf magnetic field amplitude, 2H, from equation (8), 


1 


turns out to be approximately 0.04 gauss. The optimum value of the 


peak rf current can then be calculated from the relation, 


L qT, = 2H nA (15) 
where z = | e 2nive abamperes (16) 
rf 1 
and L = the inductance of the coil in abhenries, 


n = the number of turns in the coil, 


and A = the area of the coil in on”, 


A typical value for I, can be obtained by considering one of the four 
probes mentioned in the introduction; the low frequency probe 


(4.8 to 9.1 mc/s). This coil has an inductance of 11 microhenries, 


52 turns, and area of 0.5 ae Hence 


I, =~ 1 milliampere C179 


Modulation 


As already noted, when the frequency of the rf oscillator passes 
through the resonance condition, energy is absorbed from the field of 
the coil, This absorption manifests itself by a decrease in the 
resonant impedance of the tank circuit. Thus if the tank circuit is 


fed from a constant current source, the absorption process can be observed 
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by noting a decrease in the rf voltage across it, Normally, modulation 
of the main magnetic field, Ho» is applied in order to sweep the 
resonance and obtain a periodic absorption of energy, which can then 

be easily observed on an oscilliscope. The modulation frequencies in 
practically all cases are less than 1000 cps, Hence, when the resonance 
is being swept, the rf voltage across the tank circuit is modulated 

at an audio rate, and demodulation results in the production of an 
audio signal characteristic of the absorption process taking place. 

The amplitude of the modulating field, usually supplied by a pair of 
Helmholz coils spanning the sample, should be several times the line 
width of the sample. Typically, from equation (13), the natural line 
width of a liquid sample would be of the order of two gauss (for water 
it is much less than this). Modulating fields which can be varied 

from zero to twenty gauss are commonly employed. 

The resonance condition can also be swept by using frequency 
modulation of the rf frequency. This can be done by employing a 
capacitor in the tank circuit, in addition to the tuning capacitor, 
whose capacity can be varied over small limits at an audio rate. Knoebel 
and Hahn (Kn51) constructed a frequency modulated magnetic resonance 
detector by using as the modulating capacitor a small fixed plate in 
close proximity to a sinusoidally vibrating plate, driven at 60 cps. 

A much simpler and more elegant method of producing the frequency 
modulation is realized by using in the tank circuit a solid state 


device called a varicap, which is a diode whose capacity is determined 
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by an applied bias voltage. This bias can be varied at an audio rate 
over predetermined limits to produce the desired modulation of the 
tank circuit resonant frequency. The variation in capacity AC which 


is equivalent to a value of field modulation AH is given by 
SC = 2C AH/H | (18) 


where C is the total capacity of the tuned circuit required to 
establish the resonance condition given by equation (4), w = YH.» 
For a 20 gauss sweep, AH, at a field of 2.15 kilogauss (upper limit 
of probe number one) which corresponds to value for C of 28 unf, 

AC = 0.5 wuf. The modulation circuit designed for the instrument to 


be described allows a continuous variation of AC from zero to 1.5 whrf. 


Magnetic Field Regulation 


Because the value of the magnetic field produced by an electromagnet 
is not solely dependent upon the magnet current, but also on the 
permeability of the iron (which may vary with temperature) forming the 
yoke, it is desirable to control the field by sensing any changes in 
H itself rather than by depending only on the stability of the magnet 
current power supply. This is done by utilizing the finite width of 
the proton magnetic resonance absorption line. 

If the amplitude of the modulating signal is such that the variation 
in the main magnetic field or, equivalently, the variation in the rf 


frequency in the case of frequency modulation, is greater than the line 
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width, the audio signal resulting after rf rectification has a frequency 
of 2 ow? where Wa is the modulation frequency, That is, there are two 
absorption processes for each modulation cycle, However, if the 
modulation amplitude is decreased so that only a very small part of the 
absorption curve is spanned, the resulting audio signal will have a 
frequency ow? and its amplitude will be proportional to the slope of 
the absorption curve at the point about which the field or frequency 

is varying. The maximum signal is obtained when this point is at the 
point of inflection, and the minimum, zero, when this point is at the 
maximum of the absorption curve, The phase of the audio signal 
corresponding to one side of the curve differs by 180° from that 
corresponding to the other side. 

The nature of this signal suggests the use of a phase sensitive 
detector to convert the proton signals of frequency Wa to de error 
signals, positive or negative depending on the relative phase of the 
proton signal with respect to the detector driving signal, also of 
frequency wo Such an error signal is injected into the magnet 
current power supply, thus establishing closed loop field regulation 


to within the line width of the proton sample being used. 


III, A General Description of the System 
A generalized block diagram of the system is given in figure l, 
The rf oscillator and a preamplifier form the NMR head which is located 


at the magnet whose field is to be measured and/or controlled by this 
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device. Fifty foot coaxial eaiened rei cscs AS-48) are provided to connect 
the head to the rest of the system, Remote tuning of the oscillator 
can be accomplished by means of a 0.7 rpm reversible Bodine** (Type KC1-22RC) 
motor mounted on the head and coupled to the tuning capacitor through a 
gearing system which reduces the revolution rate to approximately 0.02 rpm, 
The head itself is constructed of 1/8'"' brass plate of overall 
dimensions of 3" x 4-3/8" x 6-1/2", with separate compartments for 
the oscillator, the preamplifier, the modulation input, and the 
power line filtering components (10 millihenry" weeductors'"*** in the 
voltage lines). 
The heart of the remote, controlling part of the system is the 
twin tee oscillator. This circuit supplies the modulation voltage 
(for the varicap), the compensating signal to the preamplifier, the 
chopper drive signal, the horizontal sweep for the indicator 
oscilliscope, and a test signal to the twin tee amplifier for initial 
lining up procedures. 
The proton signal is fed directly into the Function switch, which 
then directs it to the circuits indicated in the following table, 


Table I: 


seem erence spe LCS SOL 
* Telcox Works, Greenwich, S.E. 10, England 
kk Bodine Electric Co., Chicago, Illnois, U.S.A. 

*kk Nytronics, Inc., Berkeley Heights, N.J., U.S.A. 
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The proton signal is fed directly into 
the vertical amplifier of the indicator 
oscilloscope, 


The proton signal is directed to the 
twin tee amplifier, whose output is 
fed to the oscilloscope, 


A 200 cps test signal from the twin tee 
oscillator is fed into the twin tee 
amplifier, the output of which is 
displayed on the oscilloscope, This 
position is used to test the operation 
of the twin tee amplifier, the chopper, 
and the filter circuit, 


The proton signal is fed into the twin 
tee amplifier, and the output of the 
chopper (or phase sensitive detector) 
is fed to the oscilloscope, 


When the rf frequency is sitting on the side of the absorption curve 
and the modulation amplitude is set at a small fraction of the line 
width, the proton signal has a frequency of 200 cps, that of the twin 
tee oscillator. The twin tee amplifier filters out the noise 
accompanying this proton signal and feeds it into the phase sensitive 
detector, The phase of the detector driving signal, obtained from the 
200 cps oscillator, can be made exactly in or out of phase with the 
proton signal, The detector or chopper output is then in the form of 
a full wave rectified signal. This is fed into the filter or time 


constant network where it is converted into the de error signal, Time 
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constants provided are 1, 5, 10, and 20 seconds. The 100-0-100 micro- 


ammeter provides a visual indication of the error signal. 


L¥, -<Glreuits 


A. The RF, Oscillator (Figure 2a) 

The rf oscillator employs two high frequency, silicon transistors 
(2N920, f. = 400 Mc/s) with the tank circuit in the collector circuit 
of Qa: The Varicap*, which provides the frequency modulation, is 
located in the tank circuit, It is a high Q, voltage-variable 
capacitor with a static capacitance of 10 uuf at 4 volts bias, The 
use of the transistor, Qa» in the feedback loop reduces the loading 
of the collector of Q, by the feedback elements (the 220 uuf capacitor 
and the 5K Trimpot) and, together with the 5K Trimpot, makes the 
oscillator practically independent of the emitter resistance of Qh» 
thus increasing its frequency stability and reducing the variation in 
the rf voltage level with frequency. 

The oscillator is fed from a constant current source Qo» the current 
being determined by the voltage drop across the 5.62K resistor in the 
emitter circuit of Q,> By varying the voltage of the base of Q, from 
-30 to -20 volts, the level of oscillation can be varied from 0 to 


ceoeeeneeine scatter etcetera arcane eT EEE TE EE A ILE 


* Varicap Type PC-115-10, obtained from Pacific Semiconductors, Inc., 
Hawthorne, California 
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10 volts peak-to peak over the frequency range of 4,8 to 55 Mc/s, the 
operating value being adjusted to give the optimum resonance signal, 
The peak value of the rf current (10 volts/5.62K) was chosen in 
accordance with the estimate of equation (17), allowing for a margin 
of safety. 

The above frequency range is achieved by using four different 
rf coils in the tank circuit. These will be described later. One 
setting of the feedback resistor can be found which will provide 
smooth operation of the oscillator from 4.8 to 26 mc/s, a range 
covered by three of the probes, However, when using the high frequency 
coil, that is for 26 to 55 mc/s operation, the feedback must be reduced 
slightly to obtain smooth level control over the range of this probe. 

The frequency is slightly dependent upon the value of the constant 
current driving the tank circuit so that the constant current source 
must be well stabilized and temperature effects must be reduced to a 
minimum, Hence the use of silicon transistors throughout the oscillator 
Grrcust, 

The Varicap is operated at a bias greater than or equal to 16 
volts where its capacitance vs. temperature characteristic is flat 
from =20°C to +120°C, The bias is applied through the 5,6 megohm 
resistor and is derived from the +30 volt line through a 50K panel- 
mounted potentiometer, The 200 cps modulating voltage is fed through 


the 1.5 megohm resistor. 


The choice of the low capacity varicap was made to keep to a 
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minimum the modulation of the tank circuit resonant impedance caused by 
the varying bias on this element, Such modulation produces audio 
signals which could be comparable in magnitude to the proton signal 
itself. A compensation circuit has been included in this system to 
cancel out this effect, 
The figure of merit or the Q value for this varicap at a bias 
of 16 volts and 50 mc/s is approximately 200, and exceeds 1000 at 5 Mc/s. 
The frequency stability of the oscillator was measured using a 
Hewlett Packard model 524C frequency meter, It was found that over an 
eight-hour period the oscillator was stable to one part in 10? with the 
+ 30 volts lines being supplied by transistorized power supplies. Owing 
to the sensitivity of the circuit, any disturbance on the -30 volt line 
adversely affected both the frequency stability and the signal to 
noise ratio of the proton signal. Therefore, it is highly desirable 


to use batteries for the -30 volt line of the rf oscillator circuit, 


B; The Probes 


Figure 3 shows a cross~sectional view of the probe in its container 
or shield which was designed especially for the Spectromagnetic 
Industries 50 centimeter combined broad range spectrograph and spectro- 
meter. The probe itself consists of the coil wound on a lucite form 
(which also contains the proton sample), the 1/8" brass rod, and the 
Amphenol UG-568/U connector. One end of the coil is soldered to the 
brass rod just back of the lucite form and the other to the banana 


plug in the threaded teflon end piece. The coil wire passes through 
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the end pieces which serve to keep the coil centrally psoitioned in the 
container. Teflon is used for the end piece simply because it provides 
an excellent seal to the lucite form, thus keeping the proton sample 
well contained, The banana plug fits into a blind hole drilled in the 
end of the brass container which itself serves as the return path to 
ground and as a shield for the coil, The container is permanently 
mounted in the vacuum chamber between the pole faces of the magnet with 
an O-ring vacuum seal being made at the flange. Hence, probes can be 
quickly and conveniently changed without breaking the vacuum of the 
spectrometer, The NMR head couples directly to the amphenol connector 
of the probe, of which there are four to cover the frequency range from 
4.8 to 55 mc/s as indicated in figure 3. 

In keeping with optimum coil design, the coils have Q values 
between 150 and 200, and resonant impedances of the order of 10K, 
They are 7/8" long and have a filling factor (¢ of equation (14)) of 
0.64 or, in other words, they contain 0,51 ce of the ferric nitrate 


solution mentioned earlier, 


C. The Preamplifier (Figure 2b) 

The proton resonance is detected by the 1N903 silicon diode and 
is fed through the 10K - .Qluf RC filter to a ring-of-three preamplifier 
which, with the input filter, has its upper 3db point at 1000 cps, 
Such a narrow bandwidth is in keeping with the desire to achieve the 
optimum signal to noise ratio as indicated in equation (14), The 


maximum gain of the preamplifier, given approximately by the ratio of 
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the two resistors in the emitter circuit of Q3> is 800 but this is reduced 
to about 600 by the input filter and the 1K output resistor. The noise 
at the output is approximately 5 millivolts peak-to-peak. This low 
value is achieved by having a low standing current in the first 
transistor. Thus Q, must have a high gain at low currents, a requirement 
admirably satisfied by the 2N393 transistor. 

A signal which compensates for that produced by the varying bias 
on the varicap is fed into the preamplifier through the 100K resistor 


to the emitter of Q: 
D, The Twin Tee Oscillator (Figure 4) 


The twin tee oscillator consists basically of an operational (dc 
feedback-pair) amplifier, Q, and Qo The use of such an amplifier 
ensures the de stability of the oscillator, AC negative feedback is 
applied through the twin tee notch filter for all frequencies other 
than the notch or null frequenty. Q, provides a low impedance driving 
source for the tee and Q, the high impedance load. This arrangement 
allows the values of the components in the twin tee network to be 
calculated using the standard formulas for such a filter, The loop 
gain of the circuit is essentially zero for all frequencies other than 
the notch frequency, and of course exceeds unity for that frequency, 
Slight frequency adjustments (195 - 213 cps) can be made by means of 
the 20K trimpot in the tee, and amplitude adjustments by means of the 
50K trimpot. As operating in the present system, the oscillator runs 


at its maximum output of 4 volts peak-to-peak and at 200 cps. Q, is an 
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emitter follower stage which provides isolation of the oscillator from 
the other circuits and the appropriate current to drive the modulation 
amplifier, the compensation amplifier, the horizontal amplifier, the 
chopper drive, and the twin tee amplifier when using this oscillator 
to provide a test signal for lining up procedures, 

This oscillator, unlike most other low frequency oscillators 
not employing special temperature compensating devices, is very 
stable, No detectable change in either its frequency or amplitude 


was observed when the circuit was heated to 55° Cy 


E. The Twin Tee Amplifier (Figure 5) 


The twin tee amplifier is very similar to the oscillator in 
construction. In this case, however, degeneration is provided for all 
frequencies through the 100K resistor, and additional high frequency 
degeneration through the 1000 uyf capacitor from the collector of Qa 
to ground, The tee is adjusted so that the negative feedback is 
sufficient to make the loop gain less than unity for all frequencies, The 
bandwidth is determined primarily by the setting of the 50K trimpot, and 
the notch frequency by the settings of the 20K trimpots (which should 
be the same), although they are not completely independent, The band- 
width in the present system is adjusted to 16 cps at 200 cps, where 
the peak gain is 35, The maximum input before distortion of the 
output develops is 280 millivolts peak-to-peak, Again, this amplifier 


exhibits an extremely good stability over a wide temperature range. 
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F. Modulation Amplitude Control Amplifier (Figure 6) 


This circuit is simply an operational (a dc feedback-pair) amplifier 
with a gain of 2.5, It provides a 10 volt peak-to-peak (200 cps) signal 
to the 2K helipot control on the front panel of the control unit, 

Hence, the modulation of the bias on the varicap can be varied from 


zero to + 5 volts. 


G. 200 cps Compensation Circuit (Figure 7) 


This circuit supplies up to 2 volts peak-to-peak of 200 cps 
signal of completely variable phase to the 100K resistor in the 
emitter circuit of Q, of the NMR preamplifier. Hence, any 200 cps 
signal resulting from the variation of the rf oscillator tank circuit 
impedance which is due to the varying bias on the varicap can be 
cancelled out. The boot-strapped emitter follower, Q> takes its 
input from the center tap of 20K control driven by the 4 volt peak-to-peak 
200 cps signal from the twin tee oscillator, The output of Qy is fed 
into a phase shifting network consisting of the two 1K resistors 
dividing the input to Qos the 0.5 uf capacitor feeding Q35 and the panel 
mounted potentiometer (20K) P which provides a continuous phase adjust- 
ment of 180°, Q, and Q, form a difference amplifier whose gain is 
approximately unity (given approximately by the ratio of the emitter 
resistance to the collector resistance of each transistor). The 1K 
trimpot allows the amplifier to be balanced (i.e., equality of collector 


voltages) so that maximum temperature stability can be achieved, The 
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output from the difference amplifier is taken from either collector (of 
Q or Q,) by means of a panel mounted "phase" switch and fed into the 
boot-strapped emitter follower, Qa» which provides the very high load 
impedance required by the difference amplifier. The output of Q, is 


fed into the NMR preamplifier through fifty feet of co-axial cable, 


H. Nertical Amplifier for Indicator Oscilloscope (Figure 8) 


Qy and Q, comprise a difference amplifier whose input is 
obtained from the center tap of a 10K "vertical gain" panel mounted 
control, The input to this control comes from the "Function Switch" 
on the panel, The gain of this amplifier is 13, The output from each 
collector is fed to another boot-strapped emitter follower (Q4 and Q,)» 
which in turn feeds a de fed=-back pair amplifier (Qc % and Q> Qg) 
whose gain is 15, Each de amplifier drives one vertical deflection 
plate of the oscilloscope and can supply an eighty volt peak-to-peak 
signal. Hence, an effective voltage swing of eighty volts on the 
vertical deflection plates can be obtained (vertical deflection 
sensitivity = 35 volts/inch)., The 1K trimpot in the emitter circuit 


of the difference amplifier should be adjusted to equalize the collector 


voltages of Qy and Qo. 


I, Horizontal Amplifier for Indicator Oscilloscope (Figure 9) 


This circuit is basically the same as the Vertical amplifier, H, 
In this case, however, there is a phase shifting network identical to 


that in the compensating circuit, G. The input is taken as shown from 
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the 200 cps oscillator, The 10K Horizontal Gain control and the 20K 
phase control are panel mounted, The gain of the difference 

amplifier stage is 10, and that of the de fed-back pair stage is 9, 
Each de fed=-back pair amplifier (Qe Q, and Qe» Q%) drives a horizontal 
plate of the indicator oscilloscope, The maximum effective voltage 
swing on the horizontal deflection system is approximately 90 volts 


(Horizontal deflection sensitivity = 40 volts/inch), 
J. Phase Sensitive Detector (or Chopper) Driver (Figure 10) 


This circuit consists mainly of the difference amplifier Qo» Q3 
with its input taken from the type of phase shifting network already 
described, It provides two in-phase outputs of 25 volts peak-to-peak 
of completely variable phase. These outputs are used to drive the 
chopper or phase sensitive detector circuit. The input to this 
circuit is taken directly from the twin tee oscillator (4v peak-to-peak), 
Q provides sufficient current to drive the phase shifting network, 

The output from the collector of either Q, or Q, is fed into the boot- 
strapped emitter followers Q, and Qc through the "phase reversal" 
switch mounted on the panel, Continuous phase adjustment over 180° is 


provided by the panel mounted 20K control, P. 
K. Phase Sensitive Detector (Figure 11) 


This circuit is essentially a full wave rectifier-amplifier circuit, 
Transistors er and Q, act as switches which are turned off and on by 


the Driving signal. The diodes CRl, CR2, and CR3, CR4 square off the 
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driving signal to Q and Q, respectively, CR2 and 3 limit the emitter- 
base voltage to approximately 0,7 volts, limiting the base current 
during the "on" periods to 0.7 milliamps. The biases on CR1 and CR4 
are provided by 1.3 volt mercury batteries, The 15K resistors 
preceding the diodes serve as the loads for Q, and Q. of the Driver 
circuit. 

The output of the twin tee amplifier is coupled into the 
detector circuit by means of the transformer, UTC% 0-26, Thus the 
200 cps proton signal appears at the emitters of Q, and Q 180° out 
of phase with each other. By proper adjustment of its phase, the 
driving signal can be made exactly in or out of phase with the proton 
signal, Under these conditions each transistor switch (Q) and Q,) 
half-wave rectifies the proton signal, These signals are then added 
by the adder amplifier Q. and Q%& (a dc-fed-back pair amplifier, gain = 
15) to produce an amplified, full-wave rectified proton signal at the 
output, Q, and Q,, are emitter follower isolating stages, Under no-signal 
conditions, the voltage of the bases of Q, and Q, must be adjusted for a 
minimum waveform at these points (at most two mv peak-to-peak). All the 
transistors in this circuit are dc-coupled to retain the ground reference 
for the rectified signal, The output is taken from between the 10K 
resistor and the 10K trimpot to facilitate the use of a dc-coupled, dce-— 


fed-back pair amplifier in the following time constant circuit. The 


voltage at this point is adjusted to + 15 volts, 
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L. Zime Constant Circuit (Figure 12) 


The de fed-back pair amplifier in this circuit is used to provide 
the full-wave rectified proton signal at a reference level of zero volts, 
The emitter voltage of Q, is set by the 15 volt Zener diode, and the 
emitter voltage of QQ (zero volts) is set by the 20K trimpot. The 
maximum output of this circuit to the time constant network is + 10 
volts, The 10K resistor provides the loading of Qo» and isolates the 
filter network from the amplifier, It is important to have the de 
level of Q, at zero volts so that a spurious error signal will not 
result. The filter network provides a maximum of + 5 volts at time 


constants of 1, 5, 10, or 20 seconds, 


M. The Power Input Circuit 


Figure 13 shows the circuit which provides the various voltage 
lines for the system, The input voltages, + 30 volts and - 90 volts, 
are obtained from two 100 volt, 500 milliampere, Power Design, Inc.* 
transistorized power supplies (model 105TA), The - 30 volt line for 
the oscillator section of the NMR head is taken from the battery 
circuit shown in the figure. The battery consists of two 18 volt, 
rechargeable cadmium cells in series. Zener diodes are used through- 
out to establish the various voltages. 

The indicator oscilloscope is a basic, general purpose James 
Millen** 3", rack mounted oscilloscope, model 90903, modified to allow 
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* Power Designs, Inc., 1700 Shames Drive, Westbury, N.Y. 
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signals to be applied to both plates of the vertical and the horizontal 

deflection systems, thus increasing its sensitivity by a factor of two, 
Figure 14a shows a photograph of the control system, and figure 

14b shows one of the back of the control system with its cover removed, 

Each modular circuit is lettered according to its lettering in the text. 
Figure 15a shows the NMR head coupled to one of its probes, and 


figure 15b shows the rf oscillator and preamplifier modules. 


V, Operation 


A. Lining-up Procedures 


1. Connect all the cables between the NMR head and the control 
system; that is, the power cable, the proton signal cable, the 


compensating signal cable, and the modulating signal cable, 


2. Set the Time Constant switch to the SHORT position, the rf 
oscillator amplitude control to zero, the varicap bias to half 
maximum, the modulation amplitude control to zero, and the 


compensation switch to the ‘out' position, 


3, After having set the positive and negative power supplies 
to their respective values (+ 30 volts and - 90 volts), turn on 
the power to the system and make a fine adjustment of the + 30 


volt line so that the positive and negative rails of the twin tee 


oscillator are the same, 


4, Adjust to zero the potential difference between the 
collectors in the difference amplifiers of the following 
circuits (this adjustment can be made by means of the trimpots 


in the emitter circuits of the difference amplifiers): 
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(a) the vertical amplifier (Q, and Q,) 

(b) the horizontal amplifier (Q, and Q3) 
(c) the detector driver circuit (Q,, and Q,) 
(d) the compensating circuit (Q, and Q3) 


5. Using an external oscilloscope, set the output of the 
twin tee oscillator to 4 volts peak-to-peak at 200 cps, The 20K 
trimpot provides slight frequency adjustment and the 50K trimpot 
amplitude adjustment. 


6. In the phase sensitive detector circuit, adjust to zero 
the voltages at the bases of the transistors Q3 and Q, by means of 
the 1K trimpots in the base circuits. 


7. With the Time Constant switch in the Bal position and 
observing the panel meter and/or using a voltmeter connected 
to the emitter of Q, of the Time Constant circuit, set the 
emitter voltage of this transistor to zero by adjusting the 


20K trimpot. 


8. If the zero adjustment of step 7 cannot be obtained, 
adjust the 10K trimpot in the emitter circuit of Q%& in the 


phase sensitive detector circuit and repeat steps 6 and 7, 


9, Using a 50 millivolt, 200 cps signal form the twin tee 
oscillator (set the Function switch to the No. 3 position and 
connect a jumper lead from the oscillator to the twin tee 
amplifier) check that the gain of the twin tee amplifier is 


between 35 and 40. 


10, Checking the operation of the phase sensitive detector: 
With the test signal connected as in step 9 and the time 
constant switch in the SHORT position, set the Function switch 


to the No. 4 position. This feeds the output of the chopper 
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into the indicator oscilloscope. Adjust the vertical and 
horizontal gain controls to give a display on the oscilloscope. 
Adjust the horizontal 'phase' and the chopper 'phase' to give the 
best V (or inverted) display, This indicates that the test 
signal (or the proton signal when the system is in use) is 
exactly in or out of phase with the chopper driving signal. 
This can also be checked by looking at the emitter of Q, of the 
Time Constant circuit where the waveform should be a full wave 
rectified signal (maximum peak value of + 10 volts), With the 
Function switch reset to position 3, the meter can be switched 
on to check the de output by setting the Time Constant switch 
to the Bal position, The dc output corresponding to a + 10 
volt rectified signal is + 5 volts, 


If the above steps check out, set the meter to the SHORT position, 
decrease the gain of the twin tee amplifier to zero, and remove the 
test signal jumper lead, The system is now ready for field measurement 


and control, 


Ba Magnetic Field Measurements 


1. Set the Function switch to position 1, This directs the 


signal from the NMR head to the vertical amplifier. 


2, If the approximate value of the field is known, choose 
the appropriate probe as indicated on figure 3. If probe 
No. 4 is required (i.e., for frequencies between 26 and 55 
Mc/s) refer to the discussion of the rf oscillator in section 


III for a minor adjustment of the feedback in the oscillator 


circuit. 


3. Turn up the rf amplitude until noise appears on the 


oscilloscope whose gain is set at about half maximum, 
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4, Set the modulation amplitude to half its maximum value. 
The trace on the oscilloscope will tilt because of the "spurious" 
200 cps signal being detected (see discussion of compensation 
circuit). This tilt can be cancelled out by switching in the 


compensation and adjusting its amplitude and phase. 


5. Slowly adjust the rf frequency until the proton signal 
is found. See figure 16 for the type of signal to expect. 
Center the pattern on the oscilloscope by adjusting the 
horizontal phase control, Decrease the modulation until the 
two dips are as close together as possible (any further decrease 
in the modulation will result in the disappearance of one of 
the dips). Recenter the crossover point of the two dips, and 
readjust the compensation amplitude so that the trace is 
horizontal, Under the above conditions the frequency modulation 


is sufficient to just span the absorption curve of the sample. 


6, The amplitude of the proton signal can be increased, up 
to a point, by increasing the rf amplitude. Since the rf 
amplitude affects the frequency of the oscillator, the signal 
will have to be retuned, and the compensation readjusted if the 


rf amplitude is changed from the setting of step 3. 


7. The frequency of the rf oscillator, and hence the value 


of the magnetic field, can be measured in the following way: 


Connect the signal cable from the antenna connector of the 
NMR head to a suitable preamplifier, which is to be connected 
to the frequency meter (Hewlett Packard, model 524C). After 
noting the modulation setting, decrease it to zero and measure 
the rf frequency. The value of the magnetic field can be 


obtained from the following relation: 


v = 4257,84 Ho 


where v is measured in kilocycles and Ho in gauss, 
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C, Magnetic Field Control 


1, Connect a cable from the "Error Signal" output connector 
on the control panel to the magnet power supply, with the Time 
constant switch in the BAL position, 


2. Obtain the proton signal as outlined in part B above, up 
to and including step 6, 


3. Feed the proton signal into the twin tee amplifier by 
setting the function switch to position 2, The output of this 


amplifier is now displayed on the oscilloscope. 
4, Switch out the compensation signal. 


5. Decrease the modulation level and adjust the horizontal 
phase until a circular pattern appears on the screen, Adjust 
the rf frequency slightly to optimize the amplitude of the 
circular pattern, The modulation amplitude may also have to 


be further decreased to optimize this signal, 


6. To compensate for any spurious 200 cps signal appearing 
because of the varying bias on the varicap detune the rf 
oscillator just off resonance, Increase the gain of the 
vertical amplifier, and switch in the compensation signal, 
Adjust its phase and amplitude to eliminate any circular trace 


on the screen of the oscilloscope, 


7, Retune the oscillator to the setting of step 5, When 
the optimum 200 cps signal is obtained, as indicated by the 
optimum circular pattern, the rf oscillator is sitting at the 


inflection point of the proton resonance absorption curve, 


8, To bring the detector drive signal in or out of phase 


with the proton signal adjust the detector phase control for 
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a maximum indication on the microammeter, The meter reads 


in volts and gives the magnitude of the error signal developed, 


NOTE: The meter must be in the BAL position so that hunting 


will not occur while the error signal is being obtained, 


9, With the error signal cable disconnected the phasing of 
step 8 can also be checked in the following way?! 


(a) Set the meter to SHORT position, 


(b) Set Function switch to position 4, The chopped 
signal then feeds into the vertical amplifier 
of the oscilloscope, 


(c) Adjust the horizontal phase and detector phase 
until the optimum V (or inverted) is obtained, 


(d) Set the Function switch back to position 2, and 
Time Constant switch to BAL, 


(e) Readjust horizontal phase for the circular trace 
on oscilloscope. 


(f) Reconnect the error signal cable, 


10, To start the controlling action, set the time constant to 
the appropriate time (to be determined by experiment). This 
applies the error signal to the magnet power supply. If the 
phasing between the detector drive and the proton signal is 
wrong, control will be impossible and the system will become 
regenerative. If such is the case, reverse the detector drive 


phase to establish control, 
N. Bloembergen, E. M. Purcell, and R. V. Pound, Phys. Rev., 
73, 697 (1948). 


H. W. Knoebel and E. L. Hahn, Rev. Sci. Inst., 22, 904 
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Figure 14. (a) Front view of NMR control panel 
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(b) Back view of NMR control system showing positions 
of the various circuits. 
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(a) NMR head and probe. 


Figure 15. 
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(b) RF oscillator and preamplifier modules 
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Figure 16, Typical proton resonance signals 
obtained at one kilogauss for modulation 
frequencies of (a) 200 cps, (b) 30 cps. 
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